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ABSTRACT 

We use the observed distribution of Eddington ratios as a function of supermassive black hole (BH) mass to 
constrain models of quasar/ AGN lifetimes and lightcurves. Given the observed (well constrained) AGN lumi- 
nosity function, a particular model for AGN lightcurves L{t) or, equivalently, the distribution of AGN lifetimes 
(time above a given luminosity f (> L)) translates directly and uniquely (without further assumptions) to a pre- 
dicted distribution of Eddington ratios at each BH mass. Models for self-regulated BH growth, in which feed- 
back produces a self-regulating "decay" or "blowout" phase after the AGN reaches some peak luminosity/BH 
mass and begins to expel gas and shut down accretion, make specific predictions for the lightcurves/lifetimes, 
distinct from e.g. the expected distribution if AGN simply shut down by gas starvation (without feedback) and 
very different from the prediction of simple phenomenological "light bulb" scenarios. We show that the present 
observations of the Eddington ratio distribution, spanning nearly 5 orders of magnitude in Eddington ratio, 3 
orders of magnitude in BH mass, and redshifts z = — 1, agree well with the predictions of self-regulated mod- 
els, and rule out phenomenological "light bulb" or pure exponential models, as well as gas starvation models, 
at high significance 5 a). We also compare with observations of the distribution of Eddington ratios at a 
given AGN luminosity, and find similar good agreement (but show that these observations are much less con- 
straining). We fit the functional form of the quasar lifetime distribution and provide these fits for use, and show 
how the Eddington ratio distributions place precise, tight limits on the AGN lifetimes at various luminosities, in 
agreement with model predictions. We compare with independent estimates of episodic lifetimes and use this 
to constrain the shape of the typical AGN lightcurve, and provide simple analytic fits to these for use in other 
analyses. Given these constraints, the average local BH must have gained its mass in no more than a couple 
of bright, near peak-luminosity episodes, in agreement with models of accretion triggering in interactions and 
mergers. 

Subject headings: galaxies: evolution — cosmology: theory — galaxies: active — quasars: general 



1. INTRODUCTION 

Quasars and active galactic nuclei (AGN) are among the 
most luminous, energetic, and distant objects in the Uni- 
verse. Simple integral arguments ( S oltanI 11982 ') make it 
clear that the supermassive black hole (BH) population was 
grown primarily through accretion in luminous AGN phases, 
and that the accretion luminosity released in these phases 
dominates the X-ray background and constitutes a large 
fraction of the bolometric energy prod uction of the Uni- 
verse. Comparison of e.g. the clust ering ( jCroom et al.l 120051: 
iPorciani et al.' '2004'; Hopkins et al. 2007d) and host galaxy 
properties (Bahcall et al.. .1997: CanaUzo & Stocktonl ^OOTl 



DunloD et aU l2003t iHopkins et aP l2006cl: Zakamsk aetaP 
2006t ,2008 ¥ of high and low redshift AGN and galaxies 



demonstrates that AGN are the progenitors of modern-day 
spheroids. 

Moreover, with the discovery of tight correlations be- 
tween th e masses of black holes and the veloc i ty dis - 
persion ("Ferrarese & Merritt* l2000t iGebhardt et aH l2000l) . 
masses (Magorrian et al. 1998), and perhaps m ost fundamen- 
tally binding energy or potential well depth ("Hopkins et al.' 
|20P7c,b; Aller & Richstone 2007.) of the host demonstrates 
a fundamental Unk between the growth of supermassive 
black holes and galaxy formation. A number of mod- 
els have been developed arguing that the energy or mo- 

' Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cam- 
bridge, MA 02138 

- Department of Astronomy, University of California Berkeley, Berkeley, 
CA 94720 



mentum released from an accreting supermassive black 
hole, even if only a small fraction couples to the sur- 
rounding ISM, is sufficient to halt further accretion onto 
the black hole and drive away gas, self-regulating growth 
by shutting off the quasar^ and quenching star formation 
in the galaxy and therefore allowing it to redden rapidly 
(see e .g. Ciotti & Ostriker 1997, 2001; Silk & Rees 1991 
'Burkert & Silk 2001 ; IDiMatte o et al. 2005; Hopkins et aT| 
| 2005d 2006a; Murray et al.1 ,2005; .Sazonov et al.. ,2005[ 
ISpringel et al ] l2005blah . 

One of the most basic aspects of black hole growth, and a 
powerful test of these self-regulating models for AGN evolu- 
tion, is the quasar/ AGN lifetime. Given a sufficiently well- 
known lifetime distribution (lifetime as a function of e.g. lu- 
minosity, black hole mass, and other properties), the well- 
constrained quasar luminosity function (QLF) can be empir- 
ically translated (without invoking any models or additional 
assumptions) into the triggering rate of AGN as a function of 
e.g. luminosity, BH/host galaxy/dark matter halo mass, red- 
shift, and other properties, as well as the active BH mass func- 
tion, Eddington ratio and duty cycle distributions, and differ- 
ential (mass-dependent) rate of buildup of the BH mass func- 
tion. 

Observations generally constrain q uasar lifetime s to the 
range « 10^ - 10** yr (for a review, see lMartinill2004 . These 
estimates are primarily based on demographic or integral 

^ In what follows, we use the term "quasar" somewhat loosely, as a proxy 
for high-Eddington ratio accretion activity, rather than as a reference to spe- 
cific optical properties. We use the term AGN to refer to BH accretion at all 
levels. 
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arguments which combine observations of the present-day 
population of supermassive black holes and accret i on by 
the high-redshift quasar population (e.g. . ISoltanI 



Haehnelt et al. 1998; Yu & Tremaine 2002; lYu & Lul 



1982 



2004 



Haiman et al.1 l2004; Marco ni et al.l 120041: IShankar et alJ 



20041) . or incorporate quasars into models o f galaxy eyolution 



(e.g., i Kauffmann & Haehnelt 2000; Wvith e & Loebl 1 2002 
iDi Matteo et al., ,2003i. i2004i, i2008i; Granato et al.' 70041; 
Scannapiec o & OhI 120041; iLapi et all [2006; Hopkins et al. 
2b06a/2008d; Siiacki et al. 2007) or reionization of Hell 



Sokasian et al 2002 2003; Fau cher-Giguere et alj|2008allbl; 
iMcOuinn et al.l f2009j) Resu lt s from clustering in quasar 
surveys ( e.g.. | Pprci ani et aP 120041 : |Grazian et al.' '20041; 
ICroom et al .ll2005l; iM vers et a lll2006[|20 07; Lidz et al. 2006; 
pjrciani & Norberg 2006; Sh en et aLlGOO?!; Ida Angela et alj 
1200 8; Hopkins et al. 2007d), the proximity effect in the Lya 
forest (Baitliketal. 1988; Haima n & CenI 120021; lYu & Lu 



20051: iFaucher-Giguere et al..,2008cF (but see also iLidz et al 
20071) . and the transv e rse pr o ximity effect (iJakobsen et al 
^03; Schirber et al.' 2004; Worseck & Wisotzki '2006'; 
Worseck et al. 2007; Goncalves et al. 2008) similarly suggest 
lifetimes ~ a few x 10^ yr. 

These observations, however, pertain in particular to the 
quasar lifetime - i.e. the characteristic time spent at high 
Eddington ratios/accretion rates, where much of the mass 
of a BH is accreted. Unsurpr isingly, t h e obse rvations sug- 
gest a lifetime similar to the ISalpeted (1 19641) time (the e- 
folding time for Eddington-limited black hole growth) ts = 
4.2 X 10^ (er/0.1)yr for accretion with radiative efficiency 
Cr = L/Mshc^ ~ 0.1. AGN, however, are not a homogeneous 
population, and so the lifetime is not a single number - in gen- 
eral, it should be a function of luminosity and other parame- 
ters such as BH or host mass and (possibly) various p hys- 
ical effects. As advocated by iHopkins et al.l ( 12005 dl ah . the 



AGN lifetime should properly be thought of as a luminosity- 
dependent lifetime distribution - i.e. some time f (> L) as a 
f unction of L or di fferent ial df /dlog L. 

iHopkins et al.l (l2005dllal l2006aibl) : iHopkins & Hernquist! 
(l2006ir study this luminosity-dependent lifetime/duty cycle 
distribution in hydrodynamic simulations and analytic mod- 
els of feedback-regulated BH growth, and show that such 
self-regulation leads to a generic and unique predicted form 
for the lifetime distribution. After some initial trigger that 
fuels gas inflows (such as e.g. major a nd minor mergers or 
disk instabilities; for discussion see Di Matteo et al. 2005; 
Hopkins & Hernquist 2006; Hopkins et al. 2008c, 2009a_^ 
2008al i2009cllbl: IHopkins & Hernquisti I2009t lYounger et ail 
20081) . AGN grow in approximately Eddington-limited fash- 



ion until reaching some critical mass where, if some small 
fraction of the radiant energy/momentum can couple to the 
surrounding ISM, the feedback is sufficient to halt inflows 
and/or expel gas and shut down future accretion. This "up- 
per limit" to growth is essentially an Eddington limit effective 
at the scales where the host galaxy, rather than the BH, dom- 
inates the gravitational potential, and therefore is set not by 
the details of fueling mechanisms, but instead by (relatively 
generic) global parameters of the accretion physics such as 
the host galaxy mass and AGN luminosity. 

That the resulting lifetime distribution is independent of 
fueling mechanism (i.e. is not specifically related to, for ex- 
ample, merger-induced fueling, but is a generic consequence 
of models where BH growth is self-regulated by feedback) 
has been demonstrated in the nearly identical lifetime dis- 



tributions obtained by models of fueling in major merg- 
ers, minor mergers, flyby events, bar instabilities, and ran- 
dom "stochastic" encounters with nuclear molecular clouds, 
under si milar feedback-regulated conditions (Hopkins et aT] 
2D06a b; Hopkins & Hernquist 2006t iJohansson et al.1 120091 
You ng er et a l. .2008). These various fueUng mechanisms do 
result in other important differences (in e.g. host proper- 
ties and spectral properties of observed AGN), and they will 
lead to different evolution of BHs in a cosmological sense 
(global triggering rates and their evolution varying signifi- 
cantly for the mechanisms above) and produce BHs of dif- 
ferent masses (to the extent that they produce bulges across 
a large range of masses, their self-regulating nature and the 
existence of the BH-host correlations ensures this to be the 
case ). These are discussed in more detail in other papers (see 
e.g. IHopkins & Hernquist 2009), but the important point is 
that, for a given triggering rate at some redshift and BH mass 
interval (set by some cosmological or galactic processes), 
self-regulated models predict a similar effective Eddington 
limit and lightcurve. 

Just as growth at the traditional Eddington limit leads to 
a self-similar solution for the AGN lightcurve - exponential 
growth - the expulsion of gas in this analogous limit leads to 
a self-similar lightcurve once the gas begins to be removed 
from the vicinity of the BH - a power-law decay of the form 
L cx f-(i-5-2.o) j-^j-jj^ translates into a lifetime distribu- 
tion df/dlogL with a characteristic faint-end (low-L) power- 
law like behavior: i.e. the time spent above a given luminos- 
ity /Eddington ratio (at fixed final BH mass) scales oc L^^ with 
a /3 ~ 0.6 at low-L, with a cutoff near the Eddington limit/peak 
luminosity of the system. The normalization is set by the char- 
acteristic timescales of the system - the Salpeter time and the 
(very similar) characteristic dynamical times in the central re- 
gions of the galaxy - naturally yielding a robust prediction 
of the observed quasar lifetime as well as the complete AGN 
lifetime distribution. 

This is significantly different from what is assumed in com- 
monly adopted phenomenological models, as well as other 
physical prescriptions. Quasars are often treated crudely as 
"light-bulbs" - i.e. assumed to be "on" for a fixed time (the 
"quasar lifetime") at fixed luminosity or Eddington ratio, and 
otherwise "off." In such a case, the lifetime/Eddington ra- 
tio distribution does not increase towards lower luminosities, 
but instead is strongly peaked about a characteristic high Ed- 
dington ratio/luminosity (strictly speaking, a (5-function; or 
in a Schechter-function parameterization with some scatter, 
(3 <^ 0). Similar results are obtained if one assumes that 
quasar lightcurves are pure exponentials (corresponding to 
growth at fixed Eddington ratio with either an instantaneous 
cutoff or time-mirrored exponential luminosity decay), which 
yields equal time spent per logarithmic interval in luminos- 
ity (i.e. a luminosity-independent lifetime, or /3 = Schechter 
function). 

More physically motivated but distinct models make their 
own predictions for the lifetime distribution. For example, if 
one assumes that the quasar accretion is regulated by a stan- 
dard IShakura & SunvaevI (11973 ) thin disk, and the fuel sup- 
ply is immediately removed but there is no feedback (i.e. a 
gas starvation scenario), one obtains a similarity solution for 
the accretion rate versus time that yields a lifetime distribution 
more akin to, but still significantly d istinct from that predicted 
in self-regulated models (see e.g. Y u et al.ir2005l) . 

A number of indirect tests have been proposed to break the 
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degeneracies between these models and constrain the quasar 
lifetime distribution, and the observations in these scenarios 
have thus far supported the predictions of self-regulated mod- 
els. These include: the dependence of quasar clus tering on lu- 
minosity ( Lidz et al.1 ^2006: Hopkin s et al.ll2007dl) . where ob- 
servations finding e.g. a weak dependence of clustering am- 
plitude on luminosity (at fixed redshift) support lifetime mod- 
els w i th more tirne spent at lower L jAdelberger & Steidell 
200l iCoil et alJ l2007t iMvers et al.l |lO' 



2007; da Angela et al. 
2008); the shape of the active BH mass function in various 
luminosity-selected samples (Hopkins et al. 2005bl l2006al 
|2008d), including more massive systems at lower L (in ex- 
treme cases even being peaked) rather than tracing an iden- 
tical shape to the QLF dHeckman et all 120041; iGreene & Hoi 
12007); the evolution in the faint-end QLF slope with red- 
shift (Hopkins et al. 2006b), flattening (weakly) with redshift 
jUeda et al. 2003; Hasinger et al. 2005; La Franca et al. 2005; 
Silverm an et al. 2005; Hopkins et al.. .2007 e,; , Fontanot et ah 
20071: [Silverman et al. 2008b; Sian a et al.ll2008l) as predicted 
in some self-regulated models owing to a (weak) depen- 
dence of quasar lifetime distributions on BH mass/peak lu- 
minosity; the shape o f the distribution in qu as ar host galaxy 
masses/luminosities (iHamilton et al.l l2002t iHopkins et alJ 
|2006e), similar in nature to the active BH mass function as a 
test of lifetime models; the mass functions and clustering as a 



function of mass of quasar "remnants" (i.e. bulges/spheroids, 

a re lation ; see [Ho pkins et al. 



given the observed Mbh 



2006cl 12008b 


20074 


Shankar et al. 


2009al 



Bundv et al.ll2006l;lHaimanet al.,.2007.; 
Yu&Lul 120081) 7 and the relation be- 



tween observed 

and active BH masses (Hopkins etal 



luminosity functions in different bands 
2005c; Shankar et al] 



2009bt iMerloni & Heind 120081: iBundv et al J 120081: lYu & L 
2008ir 



Although these observations are consistent with the pre- 
dictions from self-regulated models, they are indirect, and as 
such are not able to rule out alternative interpretations (in the 
case of e.g. clustering or the evolution of the faint-end QLF 
slope) or depend on additional (albeit observationally and 
physically-motivated) assumptions. Moreover, many are re- 
stricted to relatively bright Seyfert/quasar populations, where 
the models are all similar (near these luminosities, they all 
predict that the entire population must be dominated by rel- 
atively massive BHs at high Eddington ratios 0.1 — 1, and 
given the rapid growth at these Eddington ratios, the lifetime 
in this regime must be similar, comparable to the ISalpeteij 
(U964) time). Biases introduced (selecting for specific Ed- 
dington ratios, for example) will also be important in samples 
selected by e.g. broad emission lines or op tical/UV/IR col- 
ors/spectral shape (see lHopkins et al.ll2009dl) . 

The observed Eddington ratio distributions - specifically, 
the complete distribution of Eddington ratios for all BHs of 
a given mass (obscured or unobscured, luminous or under- 
luminous), however, represent a direct and powerful test of 
these models, with the predicted behavior at lower Eddington 
ratios/luminosities being strongly divergent. Given the QLF, 
a lifetime model directly and uniquely translates (without any 
additional assumptions) into a distribution of Eddington ratios 
at each BH mass, and vice versa. The requirements are de- 
manding: breaking these degeneracies necessitates large sam- 
ples, complete to all objects of a given BH mass, in a large 
volume (to constrain rare high-Eddington ratio objects), but 
sufficiently deep to measure even faint levels of AGN activity 
in those objects. Furthermore, a probe of AGN activity such 
as X-ray or narrow-line emission, robust to obscuration ef- 



fects (or the possible disappearance of the broad-emission line 
region and/or thin disk at low Eddington ratios) is important 
both to obtain a complete census of AGN activity and to avoid 
biases from e.g. luminosity or Eddington ratio-dependent ob- 
scuration/dilution/SED shapes. Fortunately, with the advent 
of wide area spectroscopic surveys such as the SDSS, it has 
become possible to constrain the Eddington ratio distribution 
at low redshifts, over a sufficiently large dynamic range to 
break the degeneracies between these models, as a function 
of various AGN and galaxy properties. 

Here, we combine a large number of observations of AGN 
Eddington ratio distributions as a function of BH and host 
galaxy mass, AGN luminosity, and redshift, in order to test 
these models and tightly constrain critical quantities such as 
the quasar lifetime as a function of luminosity, and show 
how the present observations are akeady sufficient to rule 
out, at high significance, a wide variety of alternative, simpli- 
fied physical and phenomenological AGN lifetime/lightcurve 
models. 

In § ID we compare these model predictions with observa- 
tions of the Eddington ratio distribution measured directly at 
z — Q over a range of BH mass, and inferred indirectly at 
z = 0— 1. In §[3] we similarly compare with Eddington ratio 
distributions measured not at fixed BH mass, but at fixed AGN 
luminosity, again over the range z = 0— 1. In§|4]we show 
how these observations tightly constrain physical and phe- 
nomenological models for the quasar lifetime/accretion rate 
distribution, relate to a possible dependence of quasar feed- 
back efficiency on mass, and rule out a number of alternative 
lifetime/lightcurve models. We show how the observations 
tightly constrain even general, parameterized lifetime models 
to a narrow range about the physical models, and can directly 
be converted to yield the cosmologically integrated AGN life- 
time and z = duty cycles as a function of Eddington ratio. 
In § |6] we translate these Eddington ratio/lifetime distribu- 
tion constraints to limits on the form of the "typical" AGN 
lightcurve, and discuss constraints on the "episodic" quasar 
lifetime and how, combined with the duty cycle constraints, 
this can give a bound on the number of accretion episodes per 
AGN and shape of the typical lightcurve. In §[8] we demon- 
strate the constraints from these observations on how much 
mass present-day BHs accreted in various intervals in Edding- 
ton ratio and luminosity. Finally, in §|9]we discuss our results, 
the implications of the observations for a broad range of AGN 
properties, and the prospects for future observational tests. 

For ease of comparison, we convert all observations to bolo- 
metric lum inosities given the a ppropriate bolome tric correc- 
tions fromiHopkins et alJ (l2007e.) (see also Elvis et al.lll9^4l : 
[Richards et all l2006h . We adopt a = 0.3, f^A = 0.7, 
Hq = 70kms~'Mpc~' cosmology and normalize all obser- 
vations and models appropriately (note that this generally af- 
fects only the exact normalization of quantities here, not the 
qualitative conclusions, and differences are negligible within 
the range of cosmolo gies allowed by present constraints; e.g. 
iKomatsu et ani2009h . 

2. COMPARING WITH COMPLETE EDDINGTON RATIO 
DISTRIBUTIONS 

Given the QLF ^{L\z) and some model for the quasar life- 
time or differential time at different Eddington ratios: 
df . d 



-{L\M^^) = 



t{L'>L\MsH) 



(1) 



dlogL dlogL 
it is straightforward to de-convolve and determine the Edding- 
ton ratio distribution. For example, if quasars were "light 
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bulbs" that shine either at some fixed high Eddington ratio 

L iboi 

LEdd " 1.3 X 1038ergs-i (Mbh/Mq) 



A; 



(2) 



with Aon 1 in an "on" state for a time tq and A <C 1 in an 
"off" state at other times, then the implied lifetime distribu- 
tion df/dlogL is a delta-function at Aon- The observed QLF 
in such a case is, of course, a linear translation of the active 
BH mass function (with Aon determining the re-normalization 
or shift in the horizontal axis, from Mbh to L, and the ab- 
solute value of the quasar lifetime fg determining the abun- 
dance/vertical axis). So, at a given Mbh, the Eddington ratio 
distribution will be a delta-function at Aon with a normaliza- 
tion/fractional abundance at this Aon determined by the num- 
ber of active quasars (of the corresponding luminosity) at the 
observed redshift. 

In general, in fact, so long as the quasar lifetime at a given L 
is short compared to the Hubble time (i.e. cosmological evo- 
lution in e.g. triggering rates can be ignored in the constraint), 
then the Eddington ratio distribution at a given Mbh will be 
df/dlogL (modulo a normalization reflecting the "on" frac- 
tion or tq jtu, where tu is the Hubble time at the given redshift) 
- i.e. we can trivially translate to a "duty cycle" distribution 
(fractional population at each L or A):^ 

d(5 d(5 1 df , , 

dlogA dlogL tii{z) dlogL 

By definition of the duty cycle, of course, this directly relates 
to the actual number density function of BHs at a given mass 
and luminosity or accretion rate 



$(A|Mbh) 



d«(A,MBH) 
dlogA 



: n(MBH) 



dJ 



dlogA 



(4) 



The shape of the observed Eddington ratio distribution, there- 
fore, contains information about the shape of df/dlogL in- 
dependent of either normalization, and vice versa. Such ob- 
servations hence provide a useful and direct probe of quasar 
accretion rate distributions. 

If the AGN population (e.g. triggering rates or number of 
bright objects) is relatively constant for systems of a given 
mass over the redshift range of interest (i.e. if the popula- 
tion is still growing, at least in a statistical sense, around 
the observed redshifts) and the lifetime is short relative to 
the Hubble time, then Equation |3] is applicable, and the ob- 
served Eddington ratio distribution, independent of any other 
constraints, can be directly translated to the lifetime distri- 
bution. When the triggering rate evolves strongly with red- 
shift and/or the characteristic lifetime is long compared to 
the Hubble time (e.g. for massive systems at low redshift, 
where their growth is dominated by higher-redshift periods), 
the Eddington ratio distribution at a given z is still uniquely 
predicted by a given df/dlogL model, but this must be con- 
volved over the redshift distribution of activity (i.e. given the 
observed QLF and df/dlogL, we integrate over time to de- 
termine the predicted d(5/dlogA). Even this, in practice, does 

* The translation between the distribution in logL and log A is trivial in 
terms of observations at a given BH mass; in terms of converting model pre- 
dictions between one and the other, although the relation is not completely 
trivial, we find in practice that the two are nearly equivalent, especially at 
Eddington ratios < 0.2 where most of the data with which we will compare 
lie, since some initial Mbh does not change much as a model system moves 
through a low Eddington ratio phase. We will therefore use the two inter- 
changeably in this paper as is proper for the observations, but have converted 
all physical models to the appropriate representation. 



not significantly change the direct mapping between the shape 
of the observed Eddington ratio distribution and df/dlogL; it 
mainly amounts to deriving a more correct effective "duty cy- 
cle" (normalization) than multiplying by ^ l/tn- For more 
discussion and detail s of the r elevant equations, we refer to 
.Yu & Lui (2004.) ; Hopkins et al.i (2005a, 2006 a b) ; for now, we 
note that the observed QLF $(L |z) is well-constrained at all 
the luminosities and redshifts of interest (since the relevant 
observations are primarily at low redshifts z ^ 1), and so any 
given model for df /d logL uniquely translates to an Eddington 
ratio distribution as a function of BH mass or AGN luminos- 
ity. We wish to compare these distributions to the observa- 
tions. 

In order to do so, we want to begin with a complete distri- 
bution of BH accretion rates at a given BH mass. It is impor- 
tant to do so - in models where e.g. systems of a given mass 
can have a broad range of Eddington ratios, the distribution of 
Eddington ratios at a given BH mass can be qualitatively very 
different from the distribution at a given luminosity (see §[3]l. 
Moreover, if what is desired is a distribution at all possible 
Eddington ratios, then there is no useful definition of an "ac- 
tive" AGN: what we really desire is to begin with a complete 
census of all BHs of a given mass (active or inactive at any 
level), and to measure the Eddington ratio distribution within 
this sample. 

Fortunately, the existence of a tight correlation between 
host galaxy luminosity/stellar mass/velocity dispersion and 
BH mas s make s this possible: iHeckman et al.l (120041) and 
lYu et al.l (120051) select complete samples of all SDSS galax- 
ies with a narrow range in velocity dispersion a (and cor- 
responding narrow range in BH mass), and then examine 
this sample to a limiting depth for narrow equivalent width 
A GN features (following the methodology an d classifications 
in lKauffmann et al.ll2003l; iKewlev et al.ll2006l) . Givenabolo- 
metric correction, this allows a complete census of all activity 
at a given BH mass down to some well-known limit in A. 

Moreover, the use of narrow lines is helpful for three rea- 
sons, as opposed to e.g. use of optical broad lines or an op- 
tical/IR color cut in identification of AGN and determination 
of their bolometric luminosities. First, it can probe very faint 
AGN and does not introduce much bias in terms of e.g. the 
risk of the host light in brighter systems diluting the AGN. 
Second, it allows us to include the obscured/Type 2 popula- 
tion (the_abundance of which m ay dep end on luminosi ty; see 



Bareer & Cowie 


2005; Beckmann et al. 2006; Bassani et aQ 


2006', 'Gilli et al. 


120071; IHasinsei1l2008l). Third, it should en- 



able us to identify quasars even in states of moderate radiative 
inefficiency: there is growing evidence that AGN at low Ed- 
dington ratios may transition to a radiatively inefficient state 
characterized by the absence of a thin disk. Such objects 
are still accreting but appear primarily as X-ray (because a 
hot corona survives this transition) and narrow-line sources, 
rather than broad-line or optical continuum sources (see e.g . 
|Naravan& Yi"1994'; 'Naravan etal."1995' 19961 'Meiei"2001 ; 
Maccarone et al. 2003; Yuan & Narayan 2004; Jester 20051 
iMcChntock & Remillard 2006; iCao & Xul 12007). We dis- 
cuss these distinctions in more detail in a companion pa- 
per dHopkins et al.l l2009dh . and outline how they can in- 
fluence e.g. the Eddington ratio distribution determined via 
various selection criteria. For our purposes here, however, 
the data adopted either avoid these uncertainties owing to 
their selection/identification methodology, or cover a lumi- 
nosity range where these concerns are not important (see also 
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iHeckman et ani2004l: iKauffmann et al.ll2003h . 

Figures [l]|2] show this a nalysis. We plot the Eddington ratio 
distribution determined in'Yu et al.V2005') and'Heckm an et all 
(|2b04) (see also Kauffmann & Heckman 2008), and compare 
to the theory of Hopkins et al. (2006a b, 2007e). We consider 
three versions of the estimates to show the range of theoreti- 
cal uncertainty inherent in the model, but the differences are, 
for the most part, minimal. First, we plot the iHopkins et al.l 
(|2006b) fitted dt /dlogL distribution to their typical ^ BH 
simulations, multiplied by the number density of active BHs 
at each mass (the effective duty cycle, discussed in § |4]i and 
divided by the Hubble time (since we are considering the frac- 
tion of active BHs at each A; this is shown as a the solid red 
line). This is an accurate approximation to a much more com- 
plete cosmological calculation so long as the quasar lifetime 
is short compared to the Hubble time and provided that there 
is not some strong feature in the redshift history of triggering 
(i.e. so long as the cosmological evolution around the red- 
shift of interest is relatively weak). At low A, this is simply a 
power-law; inevitably, the point where the lifetime nears the 
Hubble time will be reached, and more ful cosmological mod- 
els will turn over. 

Second, we adopt the fits from iHopkins et alJ (l2007el) to 
the observed bolometric QLF as a function of redshift, us- 
ing the lifetime distributions as a function of BH mass fitted 
in Hopkins et al. (2006b), and integrate these over redshift. 
In IHopkins et alJ (l2007e) the authors quantify the range of 
fits allowed given both the uncertainty in the lifetime model 
and the observed QLF; we show two lines that bracket this 
range (the dotted and dashed lines in Figures [T]|2]i. Again, 
these use the same model lightcurves (the power-law like 
fits to the lightcurve shape as a function of BH mass, from 
IHopkins et al.ll20()6b|) and are matched to the same QLF (the 
data compiled in IHopkins et al.ll2007el) . However, there re- 
main uncertainties in the data and degeneracies in the fit - the 
models shown bracket the range in the A distribution allowed 
by this set of models and data. The predicted distributions 
turn over at low A in both cases (albeit at slightly different 
A, reflecting these degeneracies in fitting the low-luminosity 
population), as the lifetime nears the Hubble time (obviously, 
in a full cosmological model, they must turn over, so that the 
"duty cycle" of a given BH across all Eddington ratios inte- 
grates to unity). 

The theoretical predictions all agree well, and agree re- 
markably well with the observations. We stress that these are 
predictions; no quantity has been fitted to the data points (both 
the shape and normalization of the model curves are entirely 
predicted by the papers above). 

The theoretical curves also agree with one another This is 
because, at most of the luminosities of interest, the lifetime 
is well below the Hubble time and the (especial ly low-mass) 
BHs are still growing as a population (see e.g. lHasinger et al] 
12005) . so subtleties of cosmological evolution are irrelevant 
(the only quantity of interest is df /dlogL). Moreover at low 
redshifts z < 2 the observed QLF is well-constrained, so there 
is little freedom given a df /dlogL model in what the distri- 
bution can be. The only differences appear at the highest 
masses - this is because, in a complete cosmological calcu- 
lation, most of the triggering of these systems occurred pref- 
erentially at high redshifts, so there is some pile-up at low Ed- 
dington ratios (precisely below where the condition of tg tn 
begins to be satisfied, as expected). 

Figure[2]a lso compares the Eddin gton ratio distribution de- 
termined in IHopkins et al.l (l2006dl) from the combined ob- 



served samples of iMarchesini et alj ( |2004|) and |h3 (|2002|) . 
The range in BH mass in this sample is less narrowly con- 
sti-ained, Mbh ~ 10^-^ - lO^^M©, but flie data survey ex- 
tremely faint systems in the radio and X-rays, and allow 
us to extend the observed Eddington ratio distribution from 
the already deep A ~ 10 '^ - 10"^ in the Yu et al. (2005|) 
and iHeckman et alj ( |2004|) samples by anot her two orders of 
magnitude. (Note in IHopkins et al.l (l2006dl) the authors con- 
sider possible corrections for varying radiative efficiencies 
and bolometric corrections at these luminosities, but the re- 
sulting change is for our purposes here within the error bars 
shown.) Again, the agreement is good over the entire range. 

At higher redshifts, analogous measurements are not, 
at present, available. However, we can compare with 
an alternative, albeit indirect, observational estimator It 
has been argued that the combination of radio. X-ray, 
and optical luminosities can be used to constrain the 
Eddington ratios of BHs, in a manner like that well- 
established for accreting X-ray binaries. Indeed, it is 
increasingly established that radio-loudness of AGN ap- 

fars to be a function (on average) of Eddington ratio 
arayanet al. 1995, 1996; Falcke & Biermann 1996; Meid 
2p01;'Ho 2002; Merloni et al. 2003; Marchesini et al. 200't 
Maccarone et al. 2003; Falcke et al. 2004; Merloni & Heii^ 
2007; Greene et al.i l2006i, for a review see iFender et aTl 
(2Q07J), so if this is true at higher redshift it can be 
used as at least a statistical estimator of Eddington ratio. 
iMerloni & Heinzl (12008) adopt these observations, combined 
with the measured X-ray, radio, and optical luminosity distri- 
butions of observed quasars, to constrain the Eddington ratio 
distribution as a bivariate function of BH mass and redshift. 
Their methodology allows for intrinsic scatter in these corre- 
lations, so it should be reasonably robust as long as there is 
some physical relationship between Eddington ratio and X- 
ray-radio-optical spectral shape over a wide baseline in Ed- 
dington ratio (from ^ lO^*" — 1). 

Figure [3] compares their inferred Eddington ratio distri- 
butions (A. Merloni, p rivate communication; for details see 
IMerloni & Hein3l2009h to the same predictions, at z = and 
z = 1 (at higher redshifts, IMerloni & Heini (l2008l 12009^ do 
infer Eddington ratio distributions, but the radio luminosity 
function is not directly measured and the X-ray luminosity 
function is increasingly unconstrained as well, so this re- 
lies on extrapolation of the low-redshift trends, and is not 
a direct measurement/estimate). Note that the kink in the 
Merloni & Heinz (2008) Eddington ratio distributions around 
A ~ 10^^ is sensitive to the particular assumptions about the 
form of the QLF shape and SED shape as a function of lumi- 
nosity; other attempts to i nfer this distribution hav e, however, 
seen similar features (e.g. IMarchesini et al.ir2004h . so it may 
reflect a real change in accretion properties not modeled here. 
In any case, the inferred distributions, where they overlap with 
the measurements from the observations above, agree reason- 
ably well, giving some confidence in this methodology, and 
the agreement with theoretical predictions is good. 

As discussed in § [Tj although the particular models shown 
were fitted from simulations of galaxy mergers, the results are 
comparable regardless of fueling mechanism, given a similar 
self-regulation from local AGN feedback. In fact, the pre- 
dictions and their agreement with observations in Figures [T]|3] 
rely on no information regarding fueling mechanisms - they 
simply follow from assuming a given lightcurve shape and 
matching the observed quasar luminosity function. 
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log( X = L/Le,, ) 

Fig . 1. — Distribution o f Eddington ratios at a given BH mass (at z < 0.2). We compare tlie observed distribution from SDSS naiTow-line objects lYu et alj 
|200ffl. 'Heckman et al."2004. black circles with error bars and orange diamonds, respectively) to the distribution predicted by the lightcurve/lifetime models in 
iHopkins et al. 1 2006a b, 2007e, lines). Solid red line is the simplest model prediction for a population of single triggers with Iq <^tn', black hues integrate over a 
complete cosmological history of triggering events in the model (constrained to match the observed AGN luminosity functions). Dashed and dotted lines bracket 
the model uncertainty. The lines are predictions - there are no free parameters fitted or tuned to match the observed A distributions. 



O 



O 



O 
Q. 



















= 7.9-8.1 
: (a = 176- 197) 







-8 



10 



11 















Mbh = 8.7-8.9 
_ (o = 279 -313) 




-3 -2 


-1 



log( X = L/Le,, ) 

Fig. 2. — As Figure^ continued to higher BH masses. The lower right panel shows the Eddington ratio distribution calculated in lHopkins et a f2006d) from 
the observed samples oPMarchesini et alJ 1 2004 1) andlHof (2002), with a less tightly constrained BH mass range (Mbh ~ lO'^ — 10^'^ Mq ) but extending to lower 
Eddingto n ratios. C orrections for various changes in radiative efficiency with Eddington ratio change the observed points and models within the plotted error 
bars (see lHopkins et aL,2006d) . Note that the apparent discrepancy in the models at the highest masses comes from extrapolation to unobservably small space 
densities. 



3. EDDINGTON RATIO DISTRIBUTIONS AS A FUNCTION OF 
LUMINOSITY 

At a given luminosity (as opposed to a given BH mass), the 
differences in the Eddington ratio distribution between vari- 
ous models are greatly suppressed. The reasons for this are 
obvious: at a given L and Mbh, there is only a narrow range 



of A, and given the declining number density of high-mass 
BHs, at high-L, one will increasingly be limited to the popula- 
tion of near-Eddington systems. However, over a sufficiently 
large baseline in L, differences are apparent, and (especially 
at high redshifts) luminosity-limited samples are more easily 
constructed than MBH-limited samples. We therefore consider 
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Fig. 3. — As Figure [T] but comparing th e model predictions with th e observationally inferred Eddington ratio distribution from the combination of observed 
X-ray, radio, and optical luminosities in iMerloni & Hein3 12008L|200^ and the statistical correlation between e.g. radio-lo udness and accretion rates (for a review 
see Fender et al. 2007). The results are shown over the plotted mass bins, at each labeled redshift. The sharp kinks in the lMerloni & Hem3 j2009l) distributions 
are somewhat model-dependent. 
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Fig. 4. — Predicted distribution of Eddington ratios at z ~ 1 in a narrow 
range of bolometric luminosity. Note that, because of the cut at a given lumi- 
nosity (as opposed to reflecting all A at a given Mbh as in Figure^), the dis- 
tributions are much more narrow and turn over at low Lboi . The distributions 
are reasonably approximated as log-normal, with some weakly luminosity- 
dependent skewness. 



the Eddington ratio distribution predicted as a function of lu- 
minosity at different redshifts. 

Figures |4] & H shows the results. Figure |4] demonstrates 
that at a given luminosity (unlike at a given BH mass), the 
Eddington ratio distribution is expected to be something like a 
lognormal distribution (more so if typical observational errors 
are include d), albeit with so me non-negligible skewness (as 
seen in e.g. lFine et al.ll20(j8h .^ This is because, at extremely 
low A, arbitrarily high-MeH BHs would be implied (in a bin 

' The skewness originates because the distribution of host BHs is not flat in 
mass, but increases to lower masses (following the Schechter galaxy/spheroid 
mass function). With respect to the median Eddington ratio/BH mass con- 
tributing to the observed population at a given luminosity, one therefore ex- 
pects (if Eddington ratio distributions at a given BH mass do not change 
rapidly with mass) that there will be a somewhat larger population of low- 



of fixed L), but the possible population of such systems is 
vanishing. We therefore find it convenient to approximate the 
predicted distributions as lognormal, and quantify the median 
A and 1 a dispersion (technically based on the IPV width to 
prevent bias from outliers or skewness) in the lognormal as a 
function of luminosity, at redshifts z = — 2. 

We cornpare th ese with several observational estimates. In 
iFine et afl (l2008h . the authors consider the Type 1 quasar pop- 
ulation near z w 1 in the 2dF survey, and estimate the distri- 
bution of BH masses in narrow bins of luminosity employing 
the commonly adopted virial BH mass estimators (based on 
the broad-line widths and the radius-luminosity relations in- 
ferred from reve rberation mapping of nearby AGN; see e.g. 
IVestergaard & P eterson 2006, and references therein). This 
allows them to consider the distribution of BH masses via this 
proxy as a function of luminosity down to near Seyfert lumi- 
nosities. Because the width of the distribution can be deter- 
mined without relying on the (systematically still uncertain) 
absolute normalizations of these calibrators, the authors de- 
cline to estimate absolute Eddington ratios (although a rough 
estimate su ggests they lie be tween ^0.1 — 1, as predicted 
here). Koll meier et alJ (l2006h use the same technique over a 
narrow luminosity range of Type 1 AGN in the AGES survey. 

At lower luminosities, these indicators are less useful (and 
the observations suggest the population is both more obscured 
and diluted by host galaxy light, making the virial mass esti- 
mators inaccessible). However, X-ray observations can probe 
Type 2 objects in this regime, where the optical light is dom- 
inated by the host galaxy and therefore a host galaxy stel- 
lar mass (and corresponding BH mass, adopting the observed 
Mbh — M^, relation from Marconi & Hunt 2003) can be esti- 
mated. It is well-established that in this regime the optical 
luminosity/stellar mass of the galaxy is approximately con- 
mass BHs at high Eddington ratio rather than high-mass BHs (an exponen- 
tially vanishing population) at low Eddington ratio. 
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Fig. 5. — Median Eddington ratio A and 1 a dispersion in Eddi ng ton rati os at a given narrow range in bolometric AGN luminosity (as in Figure|4] assuming 
a lognormal distribution). We compare the IHopkins et al.l j2006al lh[ I2007el) model predictions (solid lines; colors denote redshift as labeled) with the observed 
distributions inferred from the distiibution of X-ray to host luminosities in optically obscured AGN i Fiore et al. 2003; Hasinger 2008; Hickox et al. 2007, 200^, 
red diamonds) (translating host luminosity to average BH mass given the observed coiTelations), and from the distribution fitted to broad-line optical samples 
in Fine et al. (2008) from the 2dF and Kollmeier et al. (2006) from AGES, using the optical virial (line-width) BH mass estimators. We also compare with the 
predictions for pure exponential AGN lightcurves (constant Eddington ratio or exponential decay in L after some peak) and "Ught-bulb" models (where AGN 
are "on" or "o ff" with a m ass-independent narrow A distribution when on), forced to obey the (necessary) constraint of matching the observed AGN luminosity 
functions. The lHopkins et aU (,2006b) models agree well with the observations, the other models are ruled out at low-L. H i gh-L o ptical samples are not ideal for 
breaking the degeneracies between these models, although with sufficiently large dynamic range such as that in lFine et alj 120081) the distinction can be seen. 



stant while the X-ray luminosity changes, implying that at 
lower-L the X-ray luminosity function becomes increasingly 
a sequence in Eddington ratio. 

We plot the implied Eddington ratios and distribution in Ed- 
dington r atios as a function of luminosity from the sample of 
iHasingerl (^008), where the optical /?-band luminosity is con- 
verted to a stellar mass based on the age and mas s-dependent 
observed mean MjL ratios in iBell & de JongI (2000) and 
iBell et ai] ([2003 ). We have re-calculated these comparisons 
using e.g. the samples of Fiore et al. (2003) and Hickox et al. 
(l2007l l2009h and obtain the same result (various multiwave- 
length surveys have reached similar conclusions regarding 
this correlation), and find that changing the assumed host M/L 
within uncertainties makes little difference (for more discus- 
sion, see Hopkins et al. 2009d). For convenience to the com- 
parison here, we convert all the observed AGN luminosities 
to bolometric luminositi es using the bolometric corrections 
in iHopkins et al.l (l2007el) (using instead those in lElvis et alJ 
ll994tlMarconi et al.l2004tlRichards et al.l2006l makes no dif- 
ference). 

In addition to the theoretical predictions from the models 
in § |2l we contrast the results from a simple light-bulb AGN 
lifetime (in which both the mean Eddington ratio and disper- 
sion are constant) and a pure exponential model (in which 
df/dlogL is constant at A ^ 1, which, when the sample is 
cut by luminosity, does introduce some dependence on L, 



but much weaker than that predicted by the models in § |2]). 
The observations clearly prefer the steeper dependence of the 
more realistic lifetime models. Note that a large baseline in 
luminosity is needed to see the difference at high significance 
- at least ^2 — 3 orders of magnitude in L below (ideally 
more like 4 — 5 orders of magnitude). T his explains e.g. 
the we ak dependence seen in the sample of iKoUmeier et alJ 
(l2006l) . who note the weak dependence of mean Eddington 
ratio and width of the distribution on L, but have a relatively 
narrow range in L and are concentrated near and above ~ . 
Their observations are in fact entirely consistent with the 
model predictions and other observations (over larger base- 
lines) that do see such a dependence. 

4. COMPARISON WITH MODELS: WHAT DEGENERACIES ARE 
BROKEN? 

We now ask how unique these predictions are: in other 
words, can the observations distinguish between different 
models for the Eddington ratio/quasar lifetime distribution? 

Figure |6] compares the data and several simple models for 
the lifetime distribution. There are a number of commonly 
adopted forms for quasar lifetimes and lightcurves in the liter- 
ature, including the light-bulb, pure exponential growth (cor- 
responding to growth at fixed Eddington ratio and then either 
instantaneous decline or similar decay), and self-regulated de- 
cay models. The range of possibilities can be generally ap- 
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Fig. 6. — As Figure [T] but comparing the consequences of different lightcurve/lifetime models for the observed Eddington ratio distribution (generally 
parameterized by Equation|6] with power-law like slope — /3 at low A). We compare our previous model examples (red solid, with median /3 ~ 0.6) and the 
exponential (green dot-dashed, /3 = 0) and light-bulb (blue dotted, /3 <g; 0, reflecting a log-normal or 5-function distribution) models. We also show the effective 
physical upper limit (dark blue dashed, /3 = 1; /3 > 1 would result in BH growth that is divergent towards low- A). The models are all normalized to give the 
same number of objects at A > 0.1; this is required in order to match the observed quasar luminosity function where it is well-constrained (L > lO^'^ergs^'). 
Even freeing this normalization (i.e. allowing the AGN luminosity function to be significantly different from that observed), however, still does not allow the 
exponential or light-bulb models to match the observations. 



proximated by a Schechter fitting function: 

■jr~T ^'o(-r—) exp(-L/Lpeak) (5) 

dlogL Vipeak^ 

where Lpeak = V^^Edd with some 77 1. Equivalently, we can 
fit: 

^(MBH)^ro(^)^'exp(-A/,). (6) 

This allows for the fact that there must be some physi- 
cal (and observed) cutoff in the Eddington ratio distribution 
above A = 1, but permits an arbitrary power-law like behavior 
at lower Eddington ratios (which can approximate any near- 
similarity solution for self-regulated lightcurve decay, or a 
power-law-like spectrum of triggering activity, as well as a 
rapid/exponential rise fall and even sufficiently rapid rise sim- 
ilar to a light-bulb). In practice, we find that fitting an equa- 
tion of the form of Equation|6]to the data in e.g. Figures [T]|3] 
provides a useful and statistically good description of the data. 

In all mass intervals, we find a similar best-fit 77 « 0.2 — 0.4 
- the exact value reflects the particular choice of functional 
form, but the general value describes the observed and ex- 
pected cutoff in the A distribution at A ~ 1. We discuss 
variations in ry below, but find that it cannot vary widely: 
much lower 77 exponentially suppresses the number of bright 
sources, and much higher 77 implies no Eddington limit, both 
in conflict with AGN luminosity functions and direct observa- 
tions of the distribut ion of accretion rates in bright, broad-hn e 
systems (see e.g. .KoUmeier et al.ll2006t iGreene & Holl2007h . 
In what follows, we find statistically identical results fitting 
the observations with a free ij or fixed 77 = 0.4 (the degener- 
acy between e.g. 77 and /?, within the range allowed by obser- 
vations, is not strong). 

As a consequence, the shape of the distribution is primarily 
contained in the slope f3. Various lightcurve models make 



differing specific predictions for this slope: 

• Light-Bulb Models: Strictly speaking, in such a model 
df/dlogL is a delta function at the characteristic Aon; allowing 
for some finite width or measurement errors in A, this can be 
approximated as e.g. a lognormal distribution or a Schechter 
distribution with /3 ^ (i.e. a large negative jj). 

• Exponential Models: For pure exponential models (e.g. 
accretion at a constant Eddington ratio or exponential decay), 
one obtains df/dlogL ^constant, or /3 = 0. 

• Maximal low-A Accretion: At the opposite extreme, 
there is a physical limit (3 <\; because the fractional con- 
tribution to BH mass growth from a range in log(A) goes 
roughly as A x df/dlogL. For /3 > 1, the total growth is 
both formally divergent and (even if there is some cutoff at 
low-A) weighted towards the lowest-A values. Constraints 
from the observed BH mass function and thelSolt an (1982j) 
argument (see e.g. Yu & Tremaine '2002^, ISalucci et al.lll999t 
Shankar et al. 2004; Marconi et al. 2004), as well as other in- 
direct constraints dHopkins et a l. 2006d) imply that most BH 
growth cannot occur in extremely low Eddington ratio states. 

• Self-Regulated Models: In feedback-regulated models, 
the energy coupled to the ISM which halts quasar accretion 
leads to a nearly self-similar power-law like decay of the 
quasar lightcurve, L oc f~' /^. Hydrodynamic simulations of 
quasars in galaxy mergers (iHopkins et al.l2066h) suggest val- 
ues /3 ~ 0.6 for typical ~ L^: galaxie s, with a weak mass 
dependence. Analytic calculations in 'Hopkins & Hernquisll 
(5006) demonstrate that for a range of assumptions regard- 
ing the fueling mechanisms (for e.g. stochastic or secular 
ACjN fueling mechanisms), feedback coupling mechanism, 
timescale, efficiency, and gas properties around the BH, a 
range of /3 w 0.3 — 0.8 is possible; but simulations of these 
non-merger scenarios (Younger et al. 2008) suggest a similar, 
relatively narrow range of (3 independent of fueling mecha- 



10 



Hopkins & Hernquist 



nism. 
• Isolated 



Accretion Disk/Gas Starvation Models: 

lYu et all (I2005h : lYu & Lul (|2008|) show that a similar power- 
law d ecay is expected for a thin a-disk (Sh akura & SunvaevI 
Il973h abruptly cut off from any future fuel supply, when there 
is no feedback. The disk therefore slowly starves as gas ex- 
hausts by a combination of accretion and (possibly) star for- 
mation. The solution should be mass and host galaxy inde- 
pendent, with /? « 0.80 — 0.84. The exact solution depends 
(weakly) in detail on e.g. how the opacity and viscosity of the 
disk vary as a function of density and temperature, but the au- 
thors show that the entire set of solutions (corresponding to 
the allowed range for observationally reasonable a-disks) lie 
within this narrow interval in (3. In Figure |6] we omit this 
model for clarity, but the result lies between the self-regulated 
model and the physically maximal ((3 ~ I) model. 

• BHs Trace Stellar Evolution: It is possible that 
many low-luminosity systems are fueled by stellar mass 
loss from aging nuclear stellar populations. If the BH 
simply linearly traced this evolution (i.e. grew in this 
fueling-limited background with no feedback), the late-time 
accretion rate evolution would follow the stellar mass loss 
rate, giving a similar power-law solution with a steeper 
/3 pa 9- 1.0 (e.g.lN orman & Sco viUdl 19881: iLeitherer et alj 
ll999HBruzual & Charlotii2003tiCiotti & Ostrikedl2007h . 

Fitting each BH mass bin in Figures [T]|2] to an arbitrary 
function of the form in Equation|6](free to, rj, f3), we obtain the 
constraints on (3 shown in Figure]?] In fact, the constraints can 
be made stronger Any fitted A distribution, convolved with 
the BH mass function (imp Hcit in Figures THU o r taken from 
observations following e.g. iMarconi et alj|2004l) must repro- 
duce the observed AGN luminosity functions. In practice, this 
"anchors" the number density at high A > 0.1 (which dom- 
inate the QLF). Re-fitting the observed A distribution, with 
the fit constrained to also reproduce the observed QLF (taken 
h ere from the co mpilation of a large number of observations 
in lHopkins et ahi200 7e. but the results are not sensitive to the 
specific choice), we obtain the constraints in the top panel of 
Figure |7] (note that the errors shown in this case are not inde- 
pendent). 

Parametrically, if we assume f3 is independent of Mbh (or 
consider a range near '^L*), we obtain jS « 0.6 ±0.05, in good 
agreem ent with the theoretica l predictions from the simula- 
tions of iHopkins et~aLl (l2006b'). If we allow for a dependence 
of f3 on Mbh, parameterized for convenience as 

/3 = /37.5 + /?'log(MBH/10"Mo), (7) 

we obtain /Jy.s = 0.63 ± 0.04 and /?' = -0.29 ± 0.08 if we 
include the constraints from the observed QLF, and Pi 5 ^ 
0.63 ± 0.04 and /?' = -0. 11 ± 0.09 if we do not (the observed 
A distributions themselves are insufficient to determine if (3 
depends weakly on BH mass). 

Figure|8]summarizes the predictions of the different models 
above. The observations prefer a narrow range of slopes: it is 
not possible to match the observed Eddington ratio distribu- 
tion in light-bulb or exponential models (ruled out at 5 a), 
and simple stellar wind models are ruled out at 4(T. 

The isolated accretion disk model fares somewhat better: 
lYu et al.l (I2OO5I) show that it is consistent with the A distri- 
bution observed for low-mass BHs, a result we confirm here, 
but it is inconsistent with the observations of typical ~ L, and 
more massive BHs at > Aa. ,Yu & Lu (2008) show that such 
a model does improve upon the Ught-bulb model in simulta- 
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Fig. 7. — Slope /3 of the lifetime/lightcurve/Eddington ratio distribution 
fitted to tlie observations in Figures [Tj2] with a general Schechter-function 
parameterization (Equation |6). Top: Points are the maximum likelihood fit 
results for each mass bin, where the fits are constrained such that the observed 
quasar luminosity function must be reproduced. We repeat our fits using the 
quoted error bars in the observed Eddington ratio distributions (black) and 
allowing for an additional ~ 0.3 dex intrinsic unc ertainty (red). We compare 
with the model predictions from Hopkins et al. 1 2006b), determined from 
simulations allowing for a mass-dependent typical lifetime distribution or as- 
suming a constant (mass-independent) lifetime distribution (/3 0.6). Bot- 
tom: Fits repeated, but removing the constraint that the observed luminosity 
function must be reproduced (magenta star is the fit to the observations in 
the lower-right panel of Figure |2] which span a wider mass and A interval). 
The results are consistent, but less constrained. The results agree well with 
the model predictions, and if a match to the observed luminosity function is 
required, specifically favor the full mass-dependent model. In either case, an 
exponential (JS = 0) or light-bulb (/3 <C 0) model is ruled out at high signifi- 



neously reproducing the observed QLF and BH mass func- 
tion. However, as noted before, the QLF is primarily a tracer 
of high-A activity, and given that such activity dominates the 
accretion history of BHs (Soltan,1982) . these observational 
constraints alone cannot break the degeneracy between simi- 
lar, but slightly different lightcurve models L (x t~^^^ ((3 ~ 0.6 
and /3 ~ 0.8 corresponding to the self -regulated and i solated 
accretion dis k models, respectively). iHopkins et al.l (l2006al 
l2007ei l2008diib.) demonstrate that similar solutions for the 
QLF and BH mass function exist for the self -regulated model 
(unsurprisingly, the two yield very similar growth histories 
and A > 0.1 activity). As a consequence the low- A distri- 
bution represents an important means to break these degen- 
eracies, and disfavor the isolated accretion disk models. We 
do emphasize, however, that while the isolated disk model is 
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Fig. 8. — Illustration of how the observational constraints on /3 can dis- 
criminate between different physical and phenomenological quasar lifetime 
models. Blue shaded range shows the allowed slope /3 from the fits to the 
observations taken cumulatively, either assuming it is invariant (independent 
of Mbh) or allowing it to depend on Mbh as in Figure|2] Black lines com- 
pare different models for the self-regulating growth seen in hydrodynamic 
simulations (where e.g. the gas density profiles. ISM equation of state and 
phase breakdown, and prescriptions for feedback and BH accretion are var- 
ied; for details see Hopkins et al. 2006b; Hopkins & Hernquist 2006); these 
predictions agree well with the observations and are relatively independent 
of the detailed assumptions. If we include the constraint that the QLF must 
be reproduced simultaneously, the observations rule out mass-independent /3 
(corresponding to slightly more simplified analytic self-regulated models, de- 
scribed in the text). The /3 corresponding to pure exponential or phenomeno- 
logical "light-bulb" models are also shown, as is the region ruled out by the 
ISoltanI 0982) and other physical arguments. An alternative physical model, 
of a thin «-disk where the entire fuel supply is instantaneously provided from 
larger radii then cut off (i.e. large-scale gas inflows and feedback are ex- 
cluded)^ is shown (this yields fi 0.8 as the accretion disk consumes gas; 
see lYu e t al. 2005). Although a better description than the simplified phe- 
nomenological models, the observations can rule out such strictly "isolated" 
accretion disk exhaustion at > 3 cr (the constraints are even stronger, ~ 5 cr at 
high masses, if we also force the solution to match the observed QLF). Like- 
wise, assuming BH accretion simply traces stellar mass loss from evolving 
stellar populations (without feedback to regulate growth more efficiently) is 
ruled out (/3 Si 0.9 -1.0). 



ruled out at high formal significance, the absolute difference 
between it and self -regulated models is not large. That modest 
difference may reflect the inclusion of a small, but physically 
important term, for example outflows or winds from the ac- 
cretion disk that cause it to gas-exhaust slightly more steeply 
in time than otherwise. 

Feedback-regulated models agree well with observations. 
Distinguishing between feedback-regulated models, however, 
is difficult. In Figure |8] we compare several sub-classes: 
all follow a similar behav ior and reflect the basic scal- 
mgs m JsTlk&Rees d 19981) The models are derived in 
[Hopkins & Hernquist ('2006') treating feedback-driven out- 
flows as similarity solutions for expanding winds/blastwaves; 
the solutions adopt various equations of state and make differ- 
ent assumptions for the feedback coupling (whether e.g. feed- 
back instantaneously couples to the ISM or conti nues to cou- 
ple continuously throughout the event). Others ( Menci et al.l 
|2003; Granato et al. 2004) have obtained similar conclusions 
in different feedback models. We also compare the fits 
to t he hydrodynam ic simulations in iHopkins et al.l (l2006bh 
and lYounger et al.l (|2008) - the simulations give similar re- 
sults for typical spheroids, but break the strict (analytically 
assumed) self-similarity: they predict a weak mass depen- 
dence (more massive galaxies being more gas poor and bulge- 
dominated, and having somewhat more violent resonant angu- 
lar momentum transport and therefore more sharply peaked 



quasar/starburst activity as a consequence). The mass depen- 
dence amounts to a predicted /?' « —0.25 in Equation|7] The 
observations favor the self -regulated class of models, and may 
weakly favor the mass dependence predicted in simulations. 

S ome space of self - regulated models is ruled out. For Q'fja.m- 
ple. IHopkins & Hernqu ist (2006) argue that if accretion grows 
primarily out of cold instabilities "falling out of" the ex- 
panding feedback-driven blastwave (in an adiabatic gas with 
Bondi-like spherical accretion), then very steep /3 ~ 31/19 at 
early times is expected (only decaying to the /3 ~ 0.5 — 0.6 
seen in simulations in the late-time limit). Such a steep fi 
(shallow lightcurve decay) is not allowed. 

5. CONSTRAINTS ON AGN LIFETIMES 

These fits also allow us to quantify the quasar lifetime. 
There are, in fact, three different "quasar lifetimes" to which 
we could refer, which are constrained to varying degrees by 
the observations. 

(1) "Effective" Lifetimes: This is the lifetime defined in 
terms of the duty cycle - i.e. we define the effective quasar 
lifetime (at a given Eddington ratio A at a given BH mass 
A^bh) by the duty cycle d(5(A |MBH)/dlog A. Specifically, we 
invert Equation[?]to obtain 



dfeff(A|MBH) 

dlog A 



dlog A 



(8) 



or in terms of the integrated duty cycle above A (^[> A]), we 
have feff (> A) = tii{z)S{> A). This is directly determined by 
the observed Eddington ratio distributions at any redshift z, 
and is what is often referred to as the implied "AGN lifetime" 
from various observed statistics.^ 

(2) "Integrated" Lifetimes: Recall, Equation[3]is only ap- 
proximate. If one could view a BH or appropriate BH sub- 
population of final (remnant z = 0) mass Mbh over its com- 
plete history, then the true "integrated" lifetime or duty cycle 
at a given Eddington ratio or luminosity would be given by 
the appropriate integral 



,(>A) = Jsi>\\z)^^dz. 



(9) 



It is immediately clear by comparison with Equation [S] that 
the integrated lifetime is equal to the effective lifetime if the 
duty cycle d{> X\z) is relatively constant over the redshift 
range of interest. In practice, so long as (5(> A | z) has evolved 
relatively weakly between the time when most BHs of mass 
Mbh were formed and the observed redshift, one will ob- 
tain feff ~ tmi. This is seen to be the situation for relatively 
low-mass BHs corresponding to e.g. low-luminosity broad- 
line AGN (which both are observed to have relatively con- 
stant number densities and duty cycles and are still grow- 
ing, so their observed grow th corresponds to t hat during their 
"growth ep och"; see e.g . Ueda et al. 2003; Has inger et alj 
i2005; Ho pkins etal.ll2007 e; Shankai- et al. 2009b, and refer- 
ences therein). For more massive BHs, however, where S is 
significantly lower today than at high redshift, the two will 
not be the same. 

* In standard parlance, feff is defined in terms of the Hubble time tn ■ How- 
ever, strictly speaking t^ff has physical meaning (i.e. represents the time sys- 
tems spend "on" at some A at a given redshift) only if two conditions are 
met: t^ff <g; tn (the duty cycle <5 1) and host properties/triggering rates are 
varying on a timescale S> f^ff (of order tn)- Both of these are satisfied at low 
redshift, for all but the lowest-A populations, but may not be true at high red- 
shifts (imagine, for example, all hosts are "just formed" at some high redshift 
and are also all "on" - then t^ff = tn , but the real lifetime can be significantly 
shorter). 
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Nevertheless, the integrated lifetime is completely deter- 
mined by the z = Eddington ratio distributions for most rea- 
sonable physical models. The reason is that this is an integral 
constraint, and the total accretion needs to sum appropriately 
to produce a z = BH of mass Mbh- This is trivial to see 
if we temporarily consider the differential lifetime in terms 
of luminosity (rather than Eddington ratio) for a BH of final 
mass Mbh: if we have some reasonable approximation to the 
shape of the lifetime function (Equation |5]l - then the total 
Mbh must be given by the integral: 

Mbh = j Ldt = j L^^dlogL. (10) 

In detail, we should truncate the lifetime function in Equa- 
tion|5]at some minimum L where the integrated lifetime t ^tjj 
- we do so but note in practice for the typical /3 w 0.6 this 
lower L is sufficiently small that the integral has already con- 
verged. If we assume /3(Mbh) and rj are relatively constant 
(and whatever physics sets them appears, at least in simula- 
tions, to be local to the AGN and redshift independent), then 
solving this equation yields a normalization fo appropriate for 
determining the integrated quasar lifetime a bove ea ch L (for 
a more detailed derivation, see Hopkins et al.''2006ab. 

It is straightforward (although the continuity equations be- 
come somewhat more cumbersome and have to be solved nu- 
merically) to rewrite this derivation in terms of the Eddington 
ratio distribution (Equation |6]l, and the solution is similar^. 
Performing this exercise, we can approximate the integrated 
lifetimes at each logarithmic range in A by the same Equa- 
tion |6] (with 1] « 0.4), but with an "integrated" normalization 
to whose constrained value can be roughly approximated as 
1.26(1 - 1.8[/3-0.6]) X lO'^yr (near a constant - 10'^yr,i.e. 
« 2 — 3 Salpeter times, appropriate for growing a BH by an or- 
der of magnitude in mass). Integrating Equation|6]from some 
minimum A to A = 1 with this normalization, we obtain the 
required integrated lifetime in each A range. 

(3) "Episodic" Lifetimes: This is the Ufetime of an indi- 
vidual AGN/quasar "event" or some effective width in each 
peak in the quasar lightcurve. If AGN are excited to high Ed- 
dington ratio by some sort of trigger, then the time they spend 
in a given range of L or A as a consequence of just that trig- 
ger is the appropriate episodic lifetime. The duty cycle at a 
given redshift is (to lowest order) the product of the trigger- 
ing rate and episodic lifetime: for a fixed episodic lifetime, if 
the triggering rate is doubled, the effective lifetime doubles as 
well. Recall, it is the duty cycle (volume-averaged fraction) of 
BHs at a given Eddington ratio/luminosity that is observation- 
ally constrained by what we consider here. Therefore, based 
on the observations discussed thus far, we can only place an 
upper limit on the episodic lifetime (it cannot, of course, be 
longer than the integrated lifetime). We discuss the degen- 
eracies this implies and possibilities for direct constraints and 
complimentary lower limits to the episodic lifetime in § |6] 

Figure|9]compares the effective and integrated lifetimes de- 
termined from our fitting as a function of BH mass and mini- 
mum A. The effective lifetime, being essentially a duty cycle. 
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Fig. 9. — Effective quasar lifetime (black circles) above a given Edding- 
ton ratio A, i.e. the duty cycle <5(> A) tim es th e Hubble time, fitted as a 
function of BH mass to the data in Figures I1I2I The data are at low red- 
shift, so this should be thought of as the average lifetime for objects whose 
luminosity function is similar to that observed at low-z (see text). Dashed 
blue lines compare the model predictions, integrating over the cosmologi- 
cal history of triggering. The integrated quasar lifetime (red squares) deter- 
mined by integrating over the best-fit lifetime distribution to the observations 
in each bin of Mbh (applying the appropriate mass conservation equation) 
is shown for comparison. Dotted orange lines are the predicted integrated 
quasar lifetimes - reflecting the total time at each A (integrated over any 
arbitrary cosmological history) required by continuity to produce a BH at 
z = with the observed mass. The two lifetimes are similar for BHs at lower 
masses Mbh < IO'-^Mq (at which masses observations suggest the duty cy- 
cle/volume density of quasars has been relatively constant since z ^ 1^2), 
but different at the highest masses (where there has been a steep drop-off in 
quasar activity since z ~ 2; i.e. the systems did most of their growth preferen- 
tially at high redshifts and have lower average integrated lifetimes at higher 
A today). 



' In doing so, one is implicitly referring to the time spent at a given Ed- 
dington ratio near the final BH mass as observed at the appropriate redshift. 
There is no way to determine from the observations here, for example, if 
a IO'Mq BH spent many Salpeter times growing in near-Eddington lim- 
ited fashion from a small ~ IMq seed at very high redshifts, or formed 
directly as a ^ 10''Mq seed BH. More rigorously we could write this as 
t[> L/LEdd(MBH.f)]- However in terms of the Eddington ratio distribution 
that would be measured at any epoch in a relatively narrow range of BH 
mass, the two are identical. 



can be determined directly from the Eddington ratio distribu- 
tions where they cover the necessary dynamic range. The re- 
sults are as expected - roughly constant at low BH masses, 
with a steep fall in duty cycles/effective lifetimes at high- 
Mbh- The integrated lifetimes are determined from the fit- 
ted lifetime distributions (Equation |6]l, given the continuity 
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requirement. They are also as expected - for A > 0.1, for 
example, this gives an expectation that most BHs spent ~ a 
few Salpeter times at high Eddington ratio, a requirement in 
any self-consistent Eddington-limited model. They are much 
more weakly mass-dependent - since A is dimensionless and 
BH growth at fixed A exponential, we would expect that (to 
lowest order) systems of different masses should have spent 
similar time at A over the course of their evolution. At high A, 
this is true by definition; as noted above the time at a given A 
refers to time near that A near the final BH mass - if the Ed- 
dington limit is applicable, all systems must spend a similar 
time in this regime. 

Comparing the integrated and effective lifetimes yields an 
effective ratio of duty cycles: 



'int ((5(> A|z >Zobs))' 



(11) 



where {S{> A [z > Zobs)) is an appropriate weighted average S 
(weighting over the time integral in Equation|9]and by the rel- 
ative fraction of systems with mass Mbh at Zobs formed around 
each z). Roughly speaking, this yields the ratio of the duty 
cycle at the observed redshift Zobs to that during the epoch 
in which most BHs of the given mass were actively growing. 
For local low-mass BHs, this ratio is near unity; for high-mass 
BHs, is falls steeply to somewhere between ^ 1 — 10% of its 
high-redshift value, consistent with the now-standard picture 
of "downsizing" in BH growth. 

Note that the lifetimes in Figure |9] are not dire ctly analo- 
gous to those in previous works (e.g. lYu & Tremaine 2002; 
ISalucci et al.lll999HShankar et al.ll20 04: Marconi et al. 2004, 
because those assume a simplified light-bulb model, rather 
than a continuous distribution of accretion rates. 

Also, these models usually quote the total time at high Ed- 
dington ratio at all BH masses since some adopted initial con- 
dition (not the time at some Eddington ratio in a narrow range 
of BH mass, constrained directly by observations here); as 
such, the comparison requires matching the boundary condi- 
tions for seed BHs in the light-bulb models.^ Also, 5 here is 
the duty cycle above some A, and is not directly comparable 
to the "AGN fra ction" often def i ned ob servationally. 

For example. iHeckman et al.l (l2004t) point out (as seen here 
in the sa me data) that 6(> A) dec rease s strongly wit h Mbh ; 
however, ' Kauffmann et all (l2003h and iKewlev et all (l2006h 
argue from the same data that the "AGN fraction" increases 
with galaxy mass. This apparent contradiction owes to a num- 
ber of well-known trends: bulge-to-disk ratios and velocity 
dispersions drop rapidly in low-mass galaxies, so the span of 
BH mass and galaxy mass are not the same at low masses; 
specific star formation rates also increase rapidly in low-mass 
galaxies, so for the same fractional Eddington luminosity low- 
mass systems will be less AGN-dominated; and most such 
samples are AGN luminosity-limited, and so probe more mas- 
sive hosts to lower Eddington ratios. Modeling these selec- 
tion effects is non-trivial and quite sensitive to the specific 
selection effects and wavelengths of each observed sample - 
we therefore do not attempt a direct compa rison here. We 
do note, however, that the comparison of the iHeckman et al.l 

^ Specifically, we show in Figure |9] the average time spend in a given A 
range for individual BHs with a given z = mass. Some other definitions 
adopted in the literature reflect the total time for all BHs that are in or passed 
through a given BH mass interval (e.g. massive z = BHs when they were, 
earlier, at this mass). The latter lifetime definition will be higher at lower 
masses, since it includes the time local, more massive systems spent "on" at 
earlier times getting to their present-day masses. 



(12004; iK auffmann & HeckrnanI (|2008|) and lKauffmann et al.l 
(2003); K ewley et al I (l2006l) sample s indicates the consis- 
tency of these different indicators; in iHopkins et alj (l2008db 
the authors attempt to model similar selection effects to com- 
pare the self -regulated lightcurves considered here wit h AGN 
fractions in the samples observed by Kauffmann et al.l (l2003l) 
at low redshift and Erb et al. (2006); Kriek et al.l (l2007h at 
high redshift (see also iSilverman et al.ll2008a, and references 
therein), and find consistent results. 

6. TRANSLATION TO AGN LIGHTCURVES: EPISODIC VERSUS 
INTEGRATED LIFETIMES 

Given the quasar lifetime distribution, we can ask how this 
relates to some "average" quasar lightcurve L{t\MBH)- If 
the average lightcurve were monotonic, this would be triv- 
ial: inverting dr/dlogL and integrating defines L{t). For the 
Schechter function in Equation |6| this is numerically straight- 
forward but tedious; we find a statistically identical answer 
with the convenient analytic representation of the lightcurve 



A: 



1 + {t/tQ) 



1/2 



-2/13 



(12) 



where (3 is the same fitted to the lifetime distribution in Equa- 
tion |6] and 



tQ 



to- 



(13) 



/51nlO 

relative to the fo fitted in Equation |6] (the zero-point in t in 
Equation[T2]is obviously arbitrary and here is fixed to the time 
when the lightcurve is at maximum). For the median /? w 0.6, 
this gives tQ « 5.3 x 10^ yr, similar to the Salpeter time and 
characteristic dynamical times in the central regions of galax- 
ies. It is straightforward to see that the lifetime distribution 
yielded by this lightcurve is statistically a good match to the 
Schechter function fits and observations in §|2] and it is conve- 
nient for analytic models of quasar evolution. If ^ f q/2 and 
t \t\, it is trivial to treat this as time-symmetric about some 
peak, and a nearly identical lifetime distribution is obtained 
if one assumes exponential growth up to a peak luminosity 
followed by a "decay phase" given by Equation [T2l 
Equation [12] is a good match to the decline or blowout 

S)ha se of quasar evolution in feedback-regulated simulations 
Hopkins et al.l20 06b). The assumption of monotonicity need 
not be exact, provided that the (average) fractional amplitude 
of variation in the quasar luminosity is small (factor < 2 — 3) 
on timescales shorter than the integrated quasar lifetime at the 
luminosity of interest. At some point as the quasar lifetime 
approaches the Hubble time this is almost certainly not true; 
however, for the short, high-L periods this is reasonable. 

What if, however, there were two such episodes (each com- 
pressed in duration and separated by a cosmological inter- 
val >Gyr)? Figure [To] shows a few such models (arbitrarily 
re-normalized to produce the same final BH mass). First, a 
monotonic or single event case (for convenience, we assume 
the rise and fall about each "event" are symmetric). Second, 
a superposition of three such events, each with shorter dura- 
tion so that the sum df/dlogL is conserved. Third, a super- 
position of a couple of bright events and a larger number of 
lower-peak luminosity events, with a power spectrum follow- 
ing df /d log L. And fourth, a superposition of such events with 
effectively short episodic lifetimes (tg <C feft), or equivalently 
large-amplitude variability on short timescales added accord- 
ing to a power-law spectrum with slope [3. 

All of these model lightcurves reproduce the same lifetime 
distribution d//dlogL, and given df/dlogL and the observed 
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Fig. 10. — Examples of different quasar lightcurves that match the ob- 
served quasar lifetime distribution (Equation [6). The models are distin- 
guished by different episodic quasar Ufetimes and frequency of triggering. 
The integrated quasar Hfetime/duty cycle and average post-peak decay are 
the same (required to match the obs erve d A distributions) - the events decay 
after each peak following Equation [12] Different observations, for example 
probing the transverse proximity effect, are required to break the degeneracies 
between these models. Constraints from these observations at present appear 
to rule out the model in the lower right panel, but the others, with an episodic 
lifetime ~ 0.3 — 1 times the effective/integrated lifetime, are allowed. 



Eddington ratio distribution, the observed QLF will neces- 
sarily be the same. Moreover, they can all be normalized to 
give the same final BH mass (although for this final BH mass, 
they may correspond to different initial BH masses; but this 
is much smaller than the final mass and so not constrained by 
the observations considered). Measurements of the Edding- 
ton ratio distribution cannot uniquely determine the average 
lightcurve; only df/dlogL. 

In physical and observable terms, this is because the Ed- 
dington ratio distribution and corresponding df/dlogL yield 
only the effective and integrated quasar lifetimes - i.e. the 
duty cycle and total lifetime - at some L, not the episodic 
lifetime. Averaged over some redshift interval, the time t{L) 
determined from df/dlogL(L) is the total time that an object 
will be active at that L in the interval. Whether that time came 
continuously - in a single episode - or in some large number 
of shorter episodes, cannot be inferred from the observations 
considered thus far 

The episodic lifetime is of particular interest because it re- 
lates to the triggering rate of quasars. If quasars are triggered 
at a rate « and have an episodic lifetime (i.e. a duration around 
luminosity L per peak or triggering episode) of f, (~ tg if the 
event is similar to the description of Equation ITS]) , then the 
number density observed is n = nf,. If we know the observed 
number density and the effective lifetime feft, then attempting 
to infer the triggering rate we can only obtain 



" feff V ti ) ' 



(14) 



i.e. we are limited by some guess for the ratio fi/feff- The 
statement that f, = fetf is equivalent to the assumption that ev- 
ery object at a given L has had only one triggering episode 
in recent time. But if f, were smaller than feff by some ratio 
N {ti = teff/N), then the implied triggering rate for the same 
number of observed systems is higher by a factor A^. 

Assuming a physical model for what triggers quasars, we 
can use observations of quantities such as e.g. the galaxy 
merger rate to determine «, and doing so for the case of galaxy 



mergers gives the result that f, ^ feff (i-e. there are typically 
^ 1 triggers per unit Hubble time). If we wish to determine n 
without reference to a specific quasar fueling model, then we 
need the episodic lifetime f,. 

The observed rate of evolution of the quasar luminosity 
function place s an upper limit f, < 10^ yr at high luminosities 
( lMartinill2004l) . but this is already larger than the integrated 
lifetime at these luminosities. Lower limits can be d erived 
from the sizes of narrow line regions (Ben nert et alJ [2002). 
but these are short (f, > 3 x 10 yr), so are not especially 
constraining (f, ^ feff is unlikely, as it requires extremely 
high rates of triggering, on timescales faster than the rele 
vant dynamical times). Indirect measures relying on match 



ing the observed BH mass function and QLF ( Yu & Tr emaine 



, , , Ql 

'2002'; ' Salucci et al.ll 19991; IShankar et al]r2004i iMarconi et al 



2004 ) and halo occupation models which compare cluster- 
ing data dPp rciani et al. 2004; Croometal. 2005; Fin e et all 
120061; iPorciani & Norbera,2006; ,da Angela et aL,2008l) con- 
strain only the effective Ufetimes. 

One probe of episodic lifetimes is the transverse proxim- 
ity effect - if individual episodes are short, the sizes of ion- 
ized bubbles around quasars should be smaller than if indi- 
vidual episodes are long. Preliminary constraints from non- 
detections (see e.g. ISchirber et"!!] l2004t iMartinil |2004 for 
a review) suggest lifetimes f,- > 10^ yr (with a strong limit 
f, ^ lO^yr) at high luminosities (A > 0.1) where the ob- 
servations are possibl e . More recently, pote ntial detections 
jJakobsen et all 120031; iGo ncalves et al.l '2008) a nd indirect 
proxi mity effects (IWorsec k & Wisotzki 2006; Wo rseck et alJ 
l20()7h in a few objects suggest f,- « 2.5 — 5 x 10^ yr (almost 
exactly f,- = fgff ^ f™" for A > 0.1 - 0.2). Similar analy- 
sis of the Gunn-Peterson effect around high-re dshift (z > 6) 
quasa rs suggests com parable episodic lifetimes (' Baitlik et alj 
1988; Haiman & Cen 2002; Yu & Lu 2005), but is more sub- 
i ect to uncertainties in the structure of the surrounding gas 
jLidz et al. 2007). 

The lengths of relativistic jets and radio lobes also imply 
episodic lifetimes, but it is not clear that the lifetime of ra- 
dio loud activity or large-scale jet formation is the same as 
(or even correlated with) the lifetime for specific bolometri- 
cally lu minous activi ty. These observations (see e.g. Scheue^ 
1995; Blundell et al.l ll999) suggest lower-limits for f, > a few 
10*^ yr, but it is important to note that the observed systems 
with large jets in these samples are primarily at lower Edding- 
ton ratio ^0.01, so this is comparable to or a factor of a few 
smaller than the total lifetime at these luminosities. For more 
luminous FR 11 sources, a lifetime similar to that from the 
transverse proximity effect, « 2 x 10 ^ yr, has been estimated 
(Bird et al. 20 0.^) (although see also [Reynolds & BegelmanI 
ri997i:lMerloni&Heinzll2007h . 

These observations, although tentative, seem to suggest an 
episodic lifetime similar to the total lifetime, in the range 
fi/feff ~ 0.3 — 1 over reasonably high luminosity ranges A ^ 
0.01 — 1. This is consistent with simulations and cosmo- 
logical models, and similar to the characteristic timescales 
of the problem (the Salpeter time for the growth of the BH, 
4.2 X 10^ yr, and characteristic dynamical/free-fall times in 
galaxy centers, ~ 10^ — 10^ yr). Improved constraints would 
permit the extension of this comparison to lower Eddington 
ratios and enable observations to uniquely determine the trig- 
gering rates of AGN as a function of BH mass, luminosity, 
and redshift. 

7. REDSHIFT EVOLUTION: QUENCHING AND DOWNSIZING 
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Most of the observations considered here are at low red- 
shift. Some redshift evolution in Eddington ratio distributions 
must occur, as indicated in Figure |9] feff ^ fint for massive 
BHs at low redshift, so there must have been some point at 
higher redshift where the duty-cycle of high-mass BHs was 
greater (when they accreted most of their mass, around z ^ 2). 
Constraints on this evolution can be obtained from the evolu- 
tion of the quasar luminosity function - however, as discussed 
above, the QLF is primarily sensitive (especially at high red- 
shifts, where the faint end is less well-constrained) to high 
luminosities A > 0.1; likewise integral constraints from the 
shape of the BH mass function (and the requirement that the 
integrated luminosity yield the appropriate final BH mass) pri- 
marily relate to these Eddington ratios, where most mass is 
gained. 

As a consequence, the QLF constrains only the appropri- 
ate combination of /3, to, and 77 in Equation |6] such that the 
duty cycle at high Eddington ratio is reproduced. As noted 
earlier, the duty cycle is equivalent to the effective lifetime fgff 
(feff(A > 0.1)/tH{z) being the high- A duty cycle). If the Hfe- 
time/duty cycle is parameterized in the form of Equation |6l 
then it is straightforward to show that this duty cycle of inter- 
est is 

feif(A>0.1)«fo?7exp(-0.1/?7), (15) 

essentially independent of the slope parameter p (that affect- 
ing the lifetime only at lower A) and only weakly depen- 
dent on 77; for the reasonable physical range 77 ^ 0.4 — 1, feff 
changes by a factor ~ 2 — 3 (77 cannot decrease much be- 
low this with redshift, or the number density of bright ob- 
jects, observed to rise, would be exponentially suppressed). 
Various observations and synthesis models have been used 
to constrain this duty cycle - since this is just the high-A 
duty cycle, "light bulb" and more sophistica ted models yield 



nearly identical results here (see e.g. .Haehn elt et al.l 
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Yu & Tremaind 120021; lYu & Lul l2004t iHaim an et al 
Marconi et al.l 120041; IShankar et al.l 12004 )). From the fits of 



1998 



these synthesis models to the QLF, or from constructing an 
analogous simple model (assuming the lightcurves here, but 
fitting to the observed QLF from Ho pkins et al. (2007e) and 
integrating the continuity equations as in the synthesis models 
above), it can be seen that the duty cycle at fixed BH mass in- 
creases with redshift (at least from z ~ — 2) in approximate 
power-law fashion, i.e. as 



feff(A>0.1)oc(l+z)" 



(16) 



with a maximum fgff = fint ~ 10** yr for this A (obviously fgff 
cannot increase beyond this point - and indeed this reproduces 
the "flattening" in duty cycles at z > 2). 

Clearly, the duty cycle of high-mass systems must increase 
more quickly with redshift than that of low-mass systems. 
From the comparison with observations in the various syn- 
thesis models above, we note that the needed evolution can be 
crudely approximated as 



In 
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M, 



BH 



(17) 



This is not a rigorous derivation; it simply provides a use- 
ful interpolation formula that approximately reproduces the 
QLF evolution. At lowest order, this power-law increase in 
the duty cycle of high-A activity simply reflects the evolu- 
tion in the observed number density of qua sars at a given 
L — LEdd(^BH) (see e.g. Hasinger et al.ll20d5h . Second order 
corrections come from e.g. the evolving number density of 



BHs, but these depend only weakly on the lightcurve model 
in this A range. 

Given this, the simplest possible model is indeed consistent 
with the observations: a model in which quasar lightcurves ((3 
and tg, in Equation fT2ll are redshift-independent. Effectively 
this means (3 and i] in Equation |6] are fixed, as is the episodic 
lifetime of quasars - all physics local to AGN evolution are 
redshift-independent, only the triggering rate (hence the duty 
cycle and corresponding feff and fo in Equation |6]l evolve. 
This is demo nstrated in models that adopt this assumption 
dHopkins et al. 2006a, 2007e, 2008d b), but the form of the 
constraints in Equation[T6]makes it implicit. 

It is clearly important to obtain direct constraints on the Ed- 
dington ratio distributions at high redshift, particularly con- 
straints on the low-A population, in order to break the de- 
generacies between this redshift-independent model and one 
in which quasar lightcurves evolve. If there is some evo- 
lution, it may indicate a difference in fueling or feedback 
modes: one could imagine different lightcurves resulting if 
the bright, high-redshift population is fueled by violent merg- 
ers where the low-redshift population is fueled by stochas- 
tic mechanisms, bar-induced inflows, or minor mergers; or if 
feedback is not important in some populations - at high red- 
shift, feedback may act efficiently as the systems of interest 
are more bulge-dominated (leading to more sharply peaked 
lightcurves), whereas at low redshift, with activity in primar- 
ily Iow-Mbh systems with large, gas-rich disks, feedback may 
be ineffective at expelling gas content and suppressing inflows 
on large scales. 

There are some constraints on f3 and tg that can be obtained 
from low-redshift observations. As discussed above, if a pop- 
ulation is still growing (i.e. there is no sharp feature in the 
redshift history of triggering/growth) then the A distribution 
today simply reflects typical low-redshift lightcurves. How- 
ever, if there are distinct differences in the redshift history 
of triggering - if one population (at the same BH mass) is 
"quenched" (ceases growth/new triggering) at a different time 
than another - then the resulting Eddington ratio distribution 
(for the same quasar lightcurves) will not be the same at z = 0. 

Figure [TT] demonstrates this with a very simple toy model. 
Assume that AGN lightcurves are universal and redshift- 
independent, given by Equation [T2l with tg = 5 x 10^ yr and 
/3 — 0.6. Objects have luminosity L = until "triggered" 
at some time, then follow this simple lightcurve decay (f in 
Equation[T2]is the time since the trigger). The probability of a 
"trigger" as a function of time we arbitrarily parameterize as 
a Gaussian in cosmic time, rising from high redshift (z ^ 6) 
to some peak at the "quenching redshift" ZQuench and then de- 
clining to z = Specifically, 



P(ti-igger|z) oc exp| 



tH{z = 6) Lf//(zQuench) 



tH{z) 



(18) 



This is roughly chosen to correspond to the shape of the evolu- 
tion in the observed AGN luminosity density, but we empha- 
size that it is just for illustrative purposes. Figure [TT] shows 
the resulting z = distribution in A, for a set of Monte Carlo 
populations evolved in this toy model, but with a different 
quenching redshift ZQuench for each. 

Populations that quench early (high ZQuench) have all de- 
cayed for a long time, to pile up at low A and form a power- 
law distribution in A above some very low minimum A. (The 
power-law must turn over at some sufficiently low A, such that 
the duty cycle integrated over all Eddington ratios is unity; the 
specific turnover around A ^ lO^** in Figure [TT] reflects the 
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Fig. 1 1. — Top: Distribution of Eddington ratios (arb itrary units) predicted 
for populations with the same lightcurve (Equation ll2l with /3 = 0.6) but dif- 
ferent toy model triggering histories. Triggering rises with time in the same 
manner until a redshift ZQucnch, when the systems are shut down (a roughly 
Gaussian rise/fall motivated by observations). For all still-active/young sys- 
tems (low ZQiicnch): the distributions asymptote to a single distribu tion given 
by the ratio of t{> L)/tu for a single lightcurve/event (Equations ll43l with 
a turnover once t{> L) ~ tn)- For quenched/inactive systems (high ZQucnch). 
the distributions "pile up" in the long power-law tail of decayed luminosi- 
ties since the epoch of triggering. Bottom: Observed distributions from 
IKauffmann & Heckmaij jlOOSi) at fixed BH mass (~ lO'' - 10* Mq) as a 
function of stellar population age (the parameter Z)„; Z)„ < 1.5 — 1.6 cor- 
responding to star-forming/still active systems, D„ ~ 2 corresponding to "red 
and dead" systems with stellar population ages ~ tu)- The observed distribu- 
tion reflects the predicted trend: a single (redshift-independent) lightcurve is 
consistent with the observed dependence of A on age/star-forming classifica- 
tion. The "log-normal" behavior in young populations simply reflects the uni- 
versal nature of the lightcurve and inevitable turnover when ^cff{> ^) ^ tn', 
allowed to decay after quenching, this becomes the power-law-like tail in old 
populations. 



specific lightcurve and triggering functional forms assumed 
here, together with the age of the Universe. Adjusting the 
late-time lightcurve behavior or early-time triggering rates 
can shift this to lower A.) 

Populations that quench late (low ZQuench) - i-c approach 
the limit of continuous growth still occurring today - show a 
power-law behavior at high-A that directly traces the quasar 
lightcurve (as discussed above, with duty cycles simply given 
by feff ~ fint), with a turnover where fgff tn- The resulting 
distribution looks roughly log-normal, and asymptotes to the 
same distribution for all low-ZQuench (all still-active popula- 
tions). 

Note that the details in Figure [TT| are somewhat arbitrary - 
we have just chosen a toy model triggering history to high- 
light the dependence on models with strong features in that 
history; it is easy to construct others. It is possible to tune 



these histories such that, for example, the high-ZQuench pre- 
dictions continue their power-law like behavior to lower A or 
shift the "turnover" in the lognormal (low-ZQuench) regime. But 
the qualitative results are insensitive to the precise parameter- 
ization of the triggering history, provided there is a significant 
feature/shut down after a given time. 

Recently, IKauffmann & HeckmanI ( |2008|) expanded upon 
the Eddingt on ratio d istributions measured bv lHeckman et alj 
(12004) and Yu et alj (2005), and quantified the A distribution 
as a function of stellar population age (specifically the ob- 
servable parameter D„, which is a tracer of activity: D„ < 1.5 
systems being still active or < 3 Gyr old, in mean stellar pop- 
ulation age, whereas D„ ^ 2 systems are quenched with ages 
^ fn)- Figure [TT| compares the A distribution measured as a 
function of stellar population properties (for BHs with fixed 
mass 10^ — 10**Mq) at z = 0. The trends are very similar to 
those predicted (and it is not hard to imagine more detailed 
star formation history models, more precisely tuned to yield 
quantitative agreement). 

This is akeady interesting: ' Kauffmann & HeckmanI (120081) 
interpret the observed trend as an indicator of two indepen- 
dent accretion modes (a "power-law" mode, in old systems, 
and a "log-normal" mode, in young ones). This could still 
be the case, but the comparison here demonstrates that the 
observed trend is also the natural expectation of a model in 
which all AGN lightcurves are identical, but there are simply 
differences in the triggering rate corresponding to when dif- 
ferent galaxies/BH populations were "quenched" or slowed 
their growth. The uniformity of the Eddington ratio distribu- 
tion in the "log-normal" regime is not surprising in this case 
- it simply reflects the fact that quasar lightcurves are simi- 
lar and that all populations in this regime of stellar popula- 
tion age are "still active" (i.e. are continuously growing, in 
a population-averaged sense; they do not have a sharp fea- 
ture of higher-redshift activity that dominated their growth); 
the "power-law" regime is simply the "log-normal" popula- 
tion allowed to decay to lower luminosities as triggering rates 
decline with redshift. 

Moreover, these observed trends place significant con- 
straints on lightcurve evolution. Figure [12] considers the 
same model predictions, but in models where the lightcurve 
parameters in Equation [12] evolve with the "triggering red- 
shift." For convenience we parameterize the evolution sim- 
ply as (3 (X {I + z)^' and tg (x {I + z)'- , where z here refers 
to the redshift where a given trigger occurs. For /?' — —0.5 
or f' = —0.5, the predictions are similar to the no-evolution 
case, consistent with the observations. Again, we stress that 
the exact mapping to observations depends on the precise red- 
shift evolution of the triggering distribution, and more impor- 
tantly on how star formation histories in detail evolve relative 
to AGN triggering histories, needed to predict observed quan- 
tities such as D„ - however, as long as the observed qualitative 
features are intact, it is always possible to find models where 
those assumptions are adjusted to give a more precise quan- 
titative match to the observations. However, for /3' = 0.2 or 
f' = 0.5, the results are qualitatively much different - they 
do not resemble the observations. The sense of the implied 
evolution in this case would be that lightcurves become more 
shallow/extended with redshift: as a consequence, quenched 
systems would not decay sufficiently relative to observations. 
In these cases, where the features of the predicted distribu- 
tions are qualitatively different from the observations, we find 
that no tuning of the triggering rate distributions is able to 
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Fig. 12. — As Figure [TT] but in models where the median lightcurve 
(parameterized by Hfetime slope /3 and episodic Hfetime tg; Equation II 2t 
evolves with redshift. Top: Lightcurves become more shallow/extended 
with increasing redshift (/3 oc (1 +z)^'^, tg oc ( 1 + z)*^'^). It is not pos- 
sible to match the observed trends from Figure 1111 with such evolution. 
Bottom: Lightcurves become more sharply peaked with increasing redshift 
(/? oc (1 +2)""'^, tQ oc (1 +2)""''). Such evolution is consistent with the 
observed trends (and z ~ 1 observations in Figure[3j. The trend of A distribu- 
tion with age/triggering history puts significant constraints on Ughtcurve evo- 
lution: evolution towards more "quiescent" models (including stellar wind or 
isolated accretion disk modes) is disallowed; no evolution or weak evolution 
towards more violent/efficient feedback modes are consistent with observa- 
tions. 



■'fix" the disagreement with observations. 
This comparison constrains the allowed range: 



-0.5 < P', <0.05 
-0.7 < ti < 0.25 



(19) 



These are not error bars; rather, they represent the allowed 
range in whic h solutions exist that can r eprodu ce the observed 
trends from Kauff mann & HeckmanI (120081) . Within this 
range, the results are also consistent with t he observationally 
inferred Eddington ratio distributions from Merloni & Heinzl 
(12009,) at z = 1 . The constraints are non-trivial: redshift evo- 
lution must be relatively mild, and the form of evolution al- 
lowed is such that lightcurves become more sharply peaked at 
higher redshifts - the sense that might be expected if trigger- 
ing events are more violent and/or feedback is more efficient. 
Evolution in the opposite sense (evolution towards the pre- 
dictions of the stellar wind fueling or isolated accretion disk 
models discussed in § |4] or towards less efficient feedback) 
is ruled out. Together with the constraints above regarding 
evolution in the duty cycle, S = feft/fn ^ IqX dN/dt, the con- 
straints on evolution in Iq imply corresponding constraints on 
the triggering rate of independent AGN events, dN/dt. 

8. IMPLICATIONS FOR THE INTEGRATED GROWTH OF BLACK 
HOLES 

Given the observationally constrained AGN lifetime distri- 
bution. Figure [13] uses it to infer how different Eddington ra- 
tios contribute to integrated black hole growth. Essentially 
this amounts to assuming that the shape of the Eddington ratio 
distribution is well-described by Equation |6] and integrating 
over the time at each Eddington ratio to determine the frac- 
tional contribution to the final BH mass from each range in 




Fig. 13. — Contribution from different Eddington ratios to the integrated 
BH growth, given the observationally constrained typical quasar lifetime dis- 
tribution (Equation |6). Top: Differential contribution to the integrated BH 
mass from each logarithmic interval in Eddington ratio A, shown with a lin- 
ear (left) and logarithmic (right) y-axis scale. Bottom: Cumulative contribu- 
tion from all Eddington ratios less than the given A. From the observations, 
the BH growth is dominated by moderate/large Eddington ratios ~ 0.2. The 
contribution from low Eddington ratios that may be radiatively inefficient 
(A < 0.01) is small (~ 20% ), but non-negligible, in agreement with various 
independent constraints (see lHopkins et al.l2006df) . 



A. Recall from § |4] the total mass accreted must add to the 
total BH mass, and so in fractional terms the absolute nor- 
malization of the lifetimes can be factored out. Also as noted 
above, it makes little difference (because little mass growth is 
contributed by such low accretion rates) whether we truncate 
the distribution at some sufficiently low A or integrate down 
to A ^ 0. For simplicity, we consider the median best-fit life- 
time distribution with (3 « 0.6, the dependence on mass that 
may be present is sufficiently weak that the results are similar 
in the observed range. 

Figure [13] shows the fractional contribution to the final BH 
mass from each logarithmic interval in A (as well as the cu- 
mulative contribution from all Eddington ratios < A). Given 
the observations, BH growth is dominated by moderate/large 
Eddington ratios ^ 0.2. The contribution from low Edding- 
ton ratios that may be radiatively inefficient (A < 0.01) is 
small (^ 20%), but it is worth noting that it is not entirely 
negligible, and this fraction is sufficient that such populations 
could be a significant contributor to the growth in some low- 
luminosity AGN populations (for more detailed comparison, 
we refer to lHopkins et al.ll2009dl) . These expectations are in 
good agree ment with var ious independent constraints, integral 
arguments (! Soltan|[T982l) . and models for the fueling of AGN 
and buildup of the BH mass function (for a review of these 
constraints and discuss ion of the contributio n of radiatively 
inefficient sources, see iHopkins et alJl2006dL and references 
therein). 

For this reason, it is easy to construct models matching both 
the observed QLF and the z = BH mass function with a "uni- 
versal" lightcurve within the constraints devel oped here. For 
the specific lightcurve models discussed here, iHopkins et alj 
(2007e) show this explicitly (see their Figure 10, compar- 
ing the predicted BH mass function from fitting the compiled 
QLF data with the same set of light curves as those conside red 
here). Similar conclusions are reached bv lYu & Lul dlOOSl) . 

9. DISCUSSION AND CONCLUSIONS 
9.1. Overview 
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We have compared observations of the Eddington ratio dis- 
tribution as a function of BH mass, redshift, and luminosity 
with various theories. We find good agreement between the 
observations and the predictions of the self-r egulated feed- 
back models for BH gr owth and evo lution from'Hopki ns et ali 
(|2005d a,b, 2006a^; iHopkins & He rnquist (2006). The 
agreement covers the entire observed dynamic range, with ob- 
servations extending down 5—6 orders of magnitude in Ed- 
dington ratio (L/Lndd ~ 10^^ — 1), over three orders of mag- 
nitude in BH mass (Mbh ~ 10'' — IO'^Mq), corresponding to 
^ 8 orders of magnitude in luminosity, and (given indirect ob- 
servational constraints) from redshifts z '--^ — 1 . The models 
are not fitted to these observations - there are no free parame- 
ters adjustable, and the range of uncertainty between different 
versions of the models constrained to reproduce the observed 
quasar luminosity function is small (within the observational 
error bars). 

The lifetime - already as constrained in a purely empirical 
sense from the observations, is clearly not a delta function - 
i.e. quasars are not "light bulbs" - rather, it is a smooth, con- 
tinuous, and relatively steep function of luminosity /Eddington 
ratio, with more time spent at lower luminosities/Eddington 
ratios. 

We also show agreement with observed Eddington ratio dis- 
tributions as a function of AGN luminosity, but these are less 
constraining, because of the inherent Eddington ratio lim- 
its implied by the selection (wavelength-dep endent effects 
further complicate comparisons; discussed in iHopkins et al.l 
|2009d). 

9.2. Implied AGN Lifetimes and Growth Histories 

The observationally implied lifetime distribution can be 
generally parameterized (in model-independent fashion) as a 
Schechter function, with a characteristic normalization life- 
time, a turnover at Eddington ratios near ^ 1 (reflecting 
a physical limit around the Eddington limit), and a faint- 
end slope (3 (such that the luminosity-dependent lifetime 
f (> A) oc at small A). The observations favor a best-fit 
P w 0.6 ±0.05 (in agreement with feedback-regulated models, 
discussed below) for typical ~ L» galaxies and BHs. Com- 
bined with constraints from the quasar luminosity function, 
there is evidence for weak BH mass- dependence of B, as pre- 
dicted by hydrodynamic simulations (IHopkins et al.ll2006bl) . 

The observations directly yield quasar duty cycles and Ufe- 
times as a function of Eddington ratio and BH mass. In 
terms of the integrated time at a given Eddington ratio in- 
terval (around some final BH mass), the "quasar" lifetime - 
i.e. lifetime at high Eddingt o n ratio s ^0.1, is similar to that 
expected from the ISalpeterl (1 19641) time and other observa- 
tional constraints, ~ 10^ yr. However, the time at lower Ed- 
dington ratios rises rapidly - with f (A > 0.01) ^ 0.5 — 1 Gyr 
and f(A> 0.001) -l-5Gyr. 

This allows us to quantify the fractional contribution to 
present-day BH mass from various intervals in A; growth is 
dominated by large A ^ 0.2, with a small - but non-negligible 
- contribution ~ 20% from accretion at low Eddington ratios 
A < 0.01, in agreement with integral arguments (Soltan 1982) 
and various independent constraints (.Hopkins et al.. 2006d 
and references therein). 

Comparison of these integral constraints with the z = ob- 
servations (comparing "integrated" and "effective" AGN life- 
times/duty cycles) reflects the increasingly established AGN 
downsizing trend: the low-mass BH population is still grow- 
ing today, the high-mass population has shut down since its 



earlier epoch of peak activity. 

9.3. Constraints on Models of Lightcurves and Lifetimes 

There is enormous constraining power in the Eddington ra- 
tio distribution at low A in BH mass-limited samples. Already, 
the observations at z = are sufficient to limit the Edding- 
ton ratio/quasar lifetime distribution to a narrow range around 
the theoretical predictions from recent hydrodynamic simula- 
tions incorporating a self-consistent model for accretion and 
feedback. Those models predict that the self-regulated nature 
of BH accretion should lead to a relatively self-similar de- 
cay phase in AGN luminosity or Eddington ratio (regardless 
of e.g. triggering mechanisms, or the exact details of feed- 
back physics), L oc f-(i-5-2.o)^ which gives rise to the observed 
power-law-like faint behavior in the lifetime distribution with 
predicted slope f3 « 0.6. 

The observations strongly rule out (at > 5 cr significance) 
simplified models for quasar lightcurves, including: "light- 
bulb" models in which quasars turn "on" for some time at a 
fixed or relatively narrow range of accretion rates, or "expo- 
nential" models in which AGN grow at fixed Eddington ratio 
and then "shut off" rapidly. 

Models where BH accretion simply traces stellar mass loss 
are also ruled out (> 4a). Stellar mass loss may still be a 
fuel source; however, the observations argue that some pro- 
cess must further regulate AGN activity, shutting down accre- 
tion more efficiently than the simple slow starvation expected 
if BH growth directly traced stellar mass loss. 

At ^ 3 — 5 cr significance, th e observations also rule out iso- 
lated accretion disk models ( Yu et al.ll2005L and references 
therein) - i.e. accretion disks fueled rapidly but then cut off 
from a future gas supply, without feedback (such that they 
evolve by gas exhaustion). These are a considerable improve- 
ment on the light bulb model, but are still ruled out formally - 
with the sense again that some process must shut down growth 
more efficiently (possibly the addition of even mild outflows 
to the isolated disk solution) - especially for high-mass BH 
populations. 

Local observations of the Eddington ratio distribution alone 
do not strongly constrain the mass or redshift dependence of 
these properties. Combining the observations here with con- 
straints from the quasar luminosity function, observations do 
favor a weak mass-dependence in the shape of the lifetime dis- 
tribution, with the sense that f3 decreases with increasing BH 
mass (l3 - 0.6 -0.2{log[MBH/10^Mo]}). This is equivalent 
to the statement that more massive systems "shut off" more 
abruptly than low-mass systems; a trend predicted by hydro- 
dynamic models, as quasar feedback becomes relatively more 
dominant (and bulge-to-disk ratios increase while global disk 
gas fractions decrease) in more massive systems. 

Recently, Kauffmann & Heckman (2008) have shown that 
the A distribution depends on the stellar population age of the 
system. We show that the observed trends naturally arise from 
a single lightcurve of the form constrained here, as a conse- 
quence of different triggering histories. Systems which are 
still active/growing (i.e. have not "quenched" their gas supply, 
BH growth, or star formation) remain at higher accretion rates 
reflecting the median lightcurve directly, while systems which 
quenched at some earlier epoch have decayed down to the 
lower A power-law-like tail of the lightcurve distribution. The 
existence of these trends, and observ ational constraints on th e 
Eddington ratio distribution at z = 1 (iMerloni & He inz 2009). 
set some constraints on the evolution of typical lightcurves 
with redshift. 
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We find that (3 and the episodic quasar lifetime fg cannot 

evolve strongly (parameterized as /3 cx (1 and fg oc (1 + 

z)'' ■ the allowed range is -0.5 < (3'^^ < 0.05, -0.7 < t'_ < 0.25). 
Evolution towards a stellar wind or isolated accretion disk so- 
lution is strongly ruled out - rather, the sense of (mild) evo- 
lution allowed is towards more sharply peaked lightcurves at 
higher redshift, which may be expected if fueling is more vio- 
lent and/or feedback is more efficient in this regime. The evo- 
lution of duty cycles (and correspondingly "effective" AGN 
lifetimes) at the high-luminosity end (A > 0. 1) is directly con- 
strained by evolution of the quasar luminosity function - this 
limits the combination of fo and ry (but not the faint-end slope 
P), and (given some constraint on the evolution of episodic 
lifetimes fg) the triggering rate of quasars to evolve according 
to Equations [T6lfT7l increasing with redshift more rapidly in 
higher-mass BHs. 

These constraints on the AGN can be transformed into con- 
straints on the average AGN lightcurve, suggesting a charac- 
teristic power-law like decay (L oc f ^(i 5-2.o)-) sjjjyi^f j-q that 
predicted in models, after rapid growth to some peak lumi- 
nosity. The observations of the Eddington ratio distribution 
tightly limit the shape of such lightcurves, but only weakly 
constrain the characteristic timescale for a single such "event" 
(i.e. the width in time of a single "peak") - the episodic AGN 
lifetime - by setting an upper limit. 

Observational constraints from e.g. the proximity effect 
in bright quasars can independently set lower limits on the 
episodic lifetime, however, and suggest that (at least over 
the observed range for the proximity effect, namely quasar 
activity with A 0.1 — 1) the episodic lifetime is compara- 
ble to the integrated lifetime, in the range fg ~ 0.3 — l.Ofint. 
This implies that the typical massive BH has experienced no 
more than ^ a couple of bright, high Eddington ratio quasar 
episodes while near its current (z = 0) mass. This is in agree- 
ment with predictions from cosmological models that asso- 
ciate the brightest quasar activity with fueling in violent ma- 
jor mergers (Hopkins et al. 2008d), for which there are only a 
couple of events expected since high redshift. Improved con- 
straints on the episodic lifetime will allow the observed quasar 
luminosity function to be more robustly translated into the dis- 
tribution of AGN triggering rates as a function of BH and host 
mass. Extending the observational constraints on the episodic 
lifetime to lower luminosities/Eddington ratios is also impor- 
tant, as one might expect that although a typical object only 
experiences a couple of triggers to near-peak activity, it could 
have many more triggers for lower-level activity. 

9.4. Other Tests and Future Work 



Our findings agree with other (less direct) inde pendent con- 
straints. Recently, for example. I Yu & Lul (l2008h showed that 
the joint evolution of AGN luminosity functions with redshift 
favors similar lightcurves. These results (from the QLF and/or 
BH mass function and integral/continuity arguments) provide 
independent support for the constraints here, but are primar- 
ily sensitive to luminous (high- A) behavior and cannot dis- 
tinguish between similar lightcurves (with slightly different 
/3), such as the isolated accretion disk or feedback-regulated 
predictions. The lightcurve shape can also be probed by the 
dependence of AGN clustering on luminosity and shape of 
the QLF or "active" BH mass function in deep samples; the 
observations at present appear to favor feedback-regulated 
models over simplified light-bulb or exponential lightcurve 
models (see Adelberger & Steidel 2005; M vers et all 120071: 
[Greene & Ho 2007; da Angela et al. 2008), but again deeper 
observations are needed to distinguish between the proposed 
physically motivated models. 

Extending the observations of the Eddington ratio distribu- 
tion to higher redshift will greatly improve these constraints as 
well as limit possible redshift evolution in quasar light curves. 
More massive black holes will be closer to their peak growth 
at higher redshifts, allowing observations to probe the times 
of greatest interest. At z ^ 1 — 3, the A > 0. 1 end of the distri- 
bution, in broad-line luminous quasars, is already constrained 
by the ob served QLF and ap plication of the virial BH mass es- 
timators (iKoUmeier et al.|[2006t iFine et al.ll2008l) . Combined 
with these constraints, smaller volume but deep redshift sur- 
veys can be used to construct samples which are complete 
to a given BH/bulge mass, and similar narrow-line searches 
in these hosts could limit the Eddington ratio distribution at 
much lower ratios. If coverage is sufficient. X-ray data can be 
used in the same manner At z > 2, complete, large-volume 
spectroscopic or X-ray samples are not available at present. 
Here, indirect tests will remain important for the near future, 
but improved constraints on what (if any) evolution is seen at 
lower redshifts in the shape of the lifetime distributions will 
considerably inform future models and observational efforts. 
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